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Abstract 

The pp-wave/BMN limit of the AdS/CFT correspondence has exposed the Maldacena 
conjecture to a new regimen of direct tests. In one line of pursuit, finite-radius curvature 
corrections to the Penrose limit (which appear in inverse powers of the string angular mo- 
mentum J) have been found to induce a complicated system of interaction perturbations to 
string theory on the pp-wave; these have been successfully matched to corresponding cor- 
rections to the BMN dimensions of M = 4 super Yang-Mills (SYM) operators to two loops 
in the 't Hooft coupling A. This result is tempered by a well-established breakdown in the 
correspondence at three loops. Notwithstanding the third-order mismatch, we proceed with 
this line of investigation by subjecting the string and gauge theories to new and significantly 
more rigorous tests. Specifically, we extend our earlier results at 0(1/ J) in the curvature 
expansion to include string states and SYM operators with three worldsheet or i?-charge 
impurities. In accordance with the two-impurity problem, we find a perfect and intricate 
agreement between both sides of the correspondence to two-loop order in A and, once again, 
the string and gauge theory predictions fail to agree at third order. 
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1 Introduction 



Recent explorations of the AdS/CFT correspondence have focused on comparing the spec- 
trum of the type IIB superstring boosted to a near-lightcone (large angular momentum J) 
trajectory on AdS§ x S 5 with the anomalous dimensions of large -R-charge operators in M = 4 
super symmetric Yang-Mills (SYM) theory. This line of investigation was launched on the 
string side by Metsaev PjJ |2] , who showed that the string worldsheet dynamics can be re- 
duced to a system of free massive fields by taking a Penrose limit of the target geometry 
that reduces AdS$ x 5* 5 to a pp-wave. On the gauge theory side, Berenstein, Maldacena 
and Nastase [3] identified the relevant class of gauge theory operators (commonly referred 
to as BMN operators) whose dimensions can be calculated perturbatively and matched to 
the string theory energy spectrum on the pp-wave. 

Attempts to push the original results further have gone in two independent directions. In 
the gauge theory, the calculation of anomalous dimensions of BMN operators has been greatly 
simplified by Minahan and Zarembo's discovery that the problem can be mapped to that 
of computing the energies of certain integrable spin chains [3]. Based on this development, 
calculations in certain sectors of the theory have been carried out to three loops in the 
't Hooft coupling A 0EI- 1 At the same time, the quantization of the Green-Schwarz string 
in the AdS 5 x S 5 background has developed far enough to enable perturbative computations 
of the effect of worldsheet interactions on the spectrum of the string when it is boosted to 
large, but finite, angular momentum J [7J |H1 E|. These two approaches lead to different 
expansions of operator anomalous dimensions (or string eigenenergies) : on the gauge theory 
side, one naturally has an expansion in the coupling constant A which is typically exact 
in i?-charge; on the string theory side one has an expansion in inverse powers of angular 
momentum J (the dual of gauge theory i?-charge) which is exact in A. 

The expansion on the string side is difficult and has so far been carried out to 0(1/ J) for 
'two-impurity' states (i.e. states with two string oscillators excited). The resulting functions 
of the loop expansion parameter A can be compared with the large .R-charge expansion 
of two-impurity BMN operators in the gauge theory to provide new and stringent tests 
of the AdS/CFT correspondence. As mentioned above, recently developed gauge theory 
technology has made it possible to compute anomalous dimensions of certain two-impurity 
BMN operators out to three-loop order. The agreement between dual quantities is perfect out 
to two-loop order but, surprisingly, seems to break down at three loops [3 El- Exactly what 
this means for the AdS/CFT correspondence is not yet clear but, given the circumstances, 
it seems appropriate to at least look for further data on the disagreement in the hope of 
finding some instructive systematics. The subject of this paper is to pursue one possible line 
of attack in which we extend the calculations described above to higher- impurity string states 
and gauge theory operators. The extension of our two-impurity results to higher impurities 
is not a straightforward matter on either side of the correspondence and gets more complex 
as the number of impurities increases. We focus here on the three-impurity case, where we 

1 We note that the conjectural three-loop computation of [5] was solidified by field theoretic methods in 

® 
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obtain results which validate our methods for quantizing the Green-Schwarz superstring; the 
agreement with gauge theory at one and two loops is impressive, though we also confirm the 
previously observed breakdown of agreement at three-loop order. 

In section 2, we give a brief review of string quantization on AdS$ x S 5 . Since we 
use exactly the same methods and notation as in our recent treatment of the two-impurity 
problem, we can refer the reader to our long paper on that topic for a more detailed exposition 
of our methods j^j. In sections 3 and 4 we present the details of the diagonalization of 
the perturbing string worldsheet Hamiltonian on degenerate subspaces of three-impurity 
states. We give a compressed discussion of general strategy, concentrating on the aspects 
of the problem which are new to the three-impurity case. An interesting new element is 
that the non-interacting degenerate subspace breaks up into several different supersymmetry 
multiplets so that the detailed accounting of multiplicities and irrep decomposition amounts 
to a stringent test that the quantization has maintained the correct nonlinearly realized 
superconformal symmetries. Section 5 is devoted to the comparison of the string theory 
spectrum with gauge theory anomalous dimensions. Since gauge theory results have not 
been obtained elsewhere for higher-impurity operators to the order we need, we generate our 
own data by doing numerical analysis, along the lines of ^H] > of the various higher-loop spin 
chains onto which the gauge theory anomalous dimension problem has been mapped. This is 
an interesting problem in its own right and we give only a brief description of our methods, 
referring the reader to a separate methods publication jTTj for details. We are able to take 
scaling limits of the numerical spin chain analyses that allow us to make clean contact with 
the string theory results. We find perfect agreement through two-loop order and, once again, 
breakdown at three loops. Section 6 is devoted to discussion and conclusions. 

Overall, the three-impurity regime of the string theory offers a much more stringent test 
of the duality away from the full plane- wave limit. While we are unable to offer a solution to 
the disagreement with gauge theory at three loops, we can confirm that that the complicated 
interacting worldsheet theory at 0(1/ R 2 ) in the curvature expansion is properly quantized 
and correct to two loops in A. 



2 String quantization on AdS$ x S 5 : brief review 

A standard presentation of the AdS$ x S 5 metric is 

ds 2 = R 2 (-cosh 2 p dt 2 + dp 2 + sinh 2 p dttj + cos 2 6 dcj) 2 + dQ 2 + sin 2 # dfi 2 ) , (2.1) 

where R is the radius of each subspace, and dVt 2 , dVt\ are separate three-spheres. The 
spacetime radius is related to the coupling of the dual SU (N c ) gauge theory by i? 4 = A (a') 2 = 
gY M N c (a') 2 . The reparameterization 

l + ^ 2 /4 1 ~ y 2 /4 
coshp = r— cosO = —— 2.2 

simplifies the spin connection and is particularly convenient for constructing the superstring 
action in the coset formalism [EJ E] • I n these coordinates, the target-space metric takes the 
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form 



ds 2 = R 2 



±z 2 



l Z 2 



dt 2 + 



\y 2 



i + \y 2 



dzkdzk 
1 - \z 2 ) 2 



+ 



dy k 'dy k > 

1 + \v 2 ) 2 



(2.3) 



The transverse 5*0(8) breaks into 5*0(4) x 50(4), which is spanned by z 2 = z k z k with 
k = 1, . . . , 4, and y 2 



ywy k ' with k' — 5, . . . , 8. We use lightcone coordinates 
t = x + , = x + + x~/R 2 , 



leading to conjugate lightcone momenta 



- p + = A - J , -p. 



-<9 A = -- 



i? 2 



# 2 



(2.4) 



(2.5) 



where J is the conserved angular momentum conjugate to <ft and A is the conserved energy 
conjugate to t (which will eventually be identified with gauge theory operator dimensions). 
The technical benefits of this lightcone coordinate choice are explained in [HE])- When 
quantizing the string in the BMN limit, p- is held fixed while J (the angular momentum of 
the string) and R (the scale of the geometry) are both taken to be large. This is possible 
because p-R 2 = J (J2.5j) and we can regard this as either a large-J or large- R limit; we will 
pass freely between the two. 

The large- R expansion of ()2.3|) is 



ds 2 



2dx + dx~ - (x A ) 2 (dx + ) 2 + (dx 



A\2 



+ - 



1 



R 2 

+0 (l/R 4 



2y 2 dx + dx~ + -(y 4 - z 4 )(dx + f + (dx~) 2 + -z 2 dz 2 



;y 2 dy 2 



(2.6) 



When this metric is used to construct the worldsheet string action, the 0(1/ E?) terms lead 
to a quadratic (free) theory of worldsheet fields and the 0(l/i? 2 ) terms lead to quartic 
interaction terms. We will give a perturbative treatment of the effect of these interactions 
on the energy spectrum of the leading free theory. 

The Green-Schwarz action of type IIB superstring theory on AdS^ x 5 5 can be expressed as 
a non-linear sigma model constructed from supersymmetric Cartan one-forms (the index 
a is a worldsheet index here) on the coset manifold G/H = [50(4,2) x 50(6)]/[50(4, 1) x 
50(5)] [2]. The superconformal algebra of the coset manifold takes the generic form 

[B„ B v ] = f^B p [F a , B v \ = fiF p {F ai F p } = f^B, , (2.7) 



where B^ (F a ) represent bosonic (fermionic) generators. The Cartan forms and super- 
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connections L a are determined by the structure constants and f£g according to 



with L% 



TOL 



TV- 



a/3 



wah 2 (tM/2)Y 



-). 



J atl*=l 



In AdS^ x 5 5 , the Lagrangian takes the form 



(2.8) 

(2.9) 
(2.10) 



Kin 



c 



wz 



-2ie 



ab 



dtL^e^Li , 



(2.11) 
(2.12) 



where T M are 5*0(9,1) Dirac gamma matrices, tj^ is the 5*0(9,1) Minkowski metric and 
diag(l, —1). The worldsheet fermi fields 9 1 (I, J = 1,2) are two 50(9, 1) Majorana- 
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Weyl spinors of the same chirality TuO 1 = 9 1 . The K-symmetry condition YqYqO 1 = 9 1 sets 
half the components of 9 1 to zero. In constructing the Hamiltonian, it will be convenient 
to define ip = V2(d l + id 2 ): ip a is an eight-component complex spinor built out of the 16 
components of 9 1 that survive the gauge fixing. It will be useful to split components of ip 



according to their transformation under the parity operator II = 7 1 7 2 7 3 7 4 , 



where 7 a , j a are 



8 x 8 50(8) gamma matrices. 

Ultimately, we want to expand the GS worldsheet Hamiltonian in powers of 1/R, along 
the lines of (12.6)1 . in order to develop a perturbation theory of the string spectrum. To do 
this, we expand all the elements of the GS action (metric, spin connections, one-forms, etc.) 
in powers of 1/R and write the lightcone coordinates x + and x~ in terms of the physical 
transverse fields (by solving the x~ equations of motion and the covariant gauge constraints 
order-by-order). Explicit quantization rules are obtained by rewriting everything in terms 
of unconstrained canonical variables after removing the second-class constraints to which 
some of the fermionic canonical variables are subject. This awkward and rather complicated 
procedure is described in detail in [SJ |§] . The final result for the Hamiltonian density, correct 
to 0(1/ J), is 



where 
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H = H 
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(2.13) 



(2.14) 
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) 2 (x lA ) 2 + 



[X 



lA\2 



,2 



2y 



1 I i A 

x p A 



[(x A ) 2 } 2 



(2.15) 
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ff = -i{[W) + (ppO](pn^)-i(^V) 2 -i(pV) 2 + (^)(pV) 



-(^ k p)(p^ k Up') 



1 



-48 {W h 4> - pi 3k p) (pVw - P ^ k W) - u, k ^ f, k') 



(2.16) 



Hbf = — 



u . 2{ , v {VA) 2 +{x' A f + ^-z')\^'-pp') 

-\(PAX tA ){ P i>' + M ~1(pI + y' 2 - ^ 2 ) 

+ l -{z' j z k ) (^ jk ^j - frf lk p) - l -{y' r y k >) (W' fc V - P7 j ' fc 'p 
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- % -{p k Pk' - z' k y' kl )^ kk 'llp\ . 



(2.17) 



The canonical commutation relations of p A , x A and if), p mean that they can be expanded 
in bosonic creation operators a A \ fermionic creation operators b%} and their annihilation op- 
erator counterparts (where n is a mode index). Upper-case Latin letters A, B, C, ■ ■ ■ G 
1,...,8 indicate vectors of the transverse 5*0(8): these are divided into the two 50(4) 
subgroups which descend from the AdS 5 and 5 5 subspaces. Lower-case letters a, b, c, • • • G 
1, . . . , 4 label vectors in 50(4)^5, while a', b',c', - ■ ■ G 5, . . . , 8 indicate vectors in 50(4) 5 5. 
Spinors of 50(9, 1) are labeled by a, f3, 7, • • • G 1, . . . , 8. All told, there are 16 oscillators per 
mode n. The Fock space of physical string states is generated by acting on the ground state 
I J) (which carries the total angular momentum of the string state) with arbitrary combi- 
nations of the above creation operators. Matrix elements of the Hamiltonian are computed 
by expanding each of the fields and conjugate momenta of H in the mode creation and 
annihilation operators and acting on the Fock space in the obvious way. 

The parameter equivalences between the string and gauge theories are determined by 
a modified AdS/CFT dictionary which emerges in the pp-wave/BMN limit. The i?-charge 
is equated on the string side with pR 2 = J, and the large i?-charge limit corresponds to 
taking J = pJR 2 — > 00 and N c — > 00, keeping the ratio N c / J 2 fixed. The gauge theory 
coupling A = gy M N c is then replaced on the string side with the so-called modified 't Hooft 
coupling A' = g\ M N c jJ 2 . 

At leading order in J -1 , H reduces to H pp and the string theory reduces to a theory of 
eight free massive bosons and fermions. Altogether, the sixteen oscillators for each mode 
n contribute oj n = VT + n 2 A' to the energy of the string and produce a highly degenerate 
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spectrum. For example, the 256 'two- impurity' states spanned by 

a^a_ n \J) b^b_ n \J) (spacetime bosons) 

<4^\\ J ) b^a B] n \J) (spacetime fermions) (2.18) 

all have the same energy E = 1\J\ + n 2 \' in the large-J limit (note that the mode indices 
of the oscillators must sum to zero due to the level- matching constraint). This degeneracy 
is broken at the first non-leading order in J -1 by the action of the perturbing Hamiltonian 

-f^int = H^b + #ff + Hbf- 

To find the spectrum, one must diagonalize the 256 x 256 perturbation matrix schemat- 
ically represented in table H The matrix of course block-diagonalizes on spacetime bosons 



(#)int 


a^a B J n \J) 


b^X\J) 




atlb$ | J) 


( J\ a£a B n 


Hbb 


Hbf 








(J\ b a n bt n 


Hbf 


Hff 








(J\ <b« n 








Hbf 


Hbf 










Hbf 


Hbf 



Table 1: Complete Hamiltonian in the space of two-impurity string states 

and spacetime fermions, so the problem is actually 128 x 128. The matrix elements of the 
non- vanishing sub-blocks are computed by expanding each of the fields (and conjugate mo- 
menta) in Hbbi Hff an d Hbf m mode creation and annihilation operators and evaluating the 
indicated Fock space matrix elements. Finding the perturbed spectrum is then a matter of 
diagonalizing the explicit matrix constructed in this fashion. The calculation is algebraically 
tedious and requires the use of symbolic manipulation programs, but the end results are 
quite simple and verify, as previously described, a perfect match of the string spectrum to 
the gauge theory results out to two-loop order. 

Our purpose in this paper is to work out the analogous results for higher-impurity states 
on both sides of the correspondence. In fact, we will limit our detailed considerations to the 
three-impurity problem, since it presents many interesting complications as compared to the 
two-impurity case, and it is not clear that any useful illumination will come from studying 
yet higher-impurity states. 

3 Three-impurity string spectrum: leading order in A' 

The three-impurity Fock space block-diagonalizes into separate spacetime fermion and space- 
time boson sectors. The bosonic sector contains states that are purely bosonic (composed of 
three bosonic string oscillators) and states with bi-fermionic components: 

afa B ^a^\J) of b? fof | J) . (3.1) 
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Pure boson states are mixed by the bosonic sector of the Hamiltonian Hbb, while states 
with bi-fermionic excitations are mixed both by the purely fermionic Hamiltonian Hpp and 
the bose-fermi sector Hbf- The sector of spacetime fermion states is composed of purely 
fermionic excitations and mixed states containing two bosonic oscillators: 

afWlJ) a?a?bf\J) . (3.2) 

Pure fermion states are acted on by Hff, and mixed states with bosonic excitations are acted 
on by Hbb and Hbf- This block diagonalization of the perturbing Hamiltonian is displayed 
schematically in table 121 



(H)i n t 


a A ^a B ^a c ^\J) 




6 Q tfr/3tfr7t |J) 


a A] a B] h a\\j^ 


(J\ a A a B a c 


Hbb 


Hbf 








(J\ a A b a bP 


Hbf 


Hff + Hbf 








(J\ b a tfW 








Hff 


Hbf 


(J\ a A a B b a 








Hbf 


Hbb + Hbf 



Table 2: Three-impurity string states 

The three-impurity string states are subject to the usual level-matching condition on the 
mode indices: q + r + s = 0. There are two generically different solutions of this constraint: 
all mode indices different (q ^ r ^ s) and two indices equal (eg. q = r = n, s = — 2n). In the 
inequivalent index case, there are 16 3 = 4, 096 degenerate states arising from different choices 
of spacetime labels on the mode creation operators. In the case of two equivalent indices, the 
dimension of the degenerate subspace is half as large (there are fewer permutations on mode 
indices that generate linearly independent states). The two types of basis break up into 
irreducible representations of PSU(2, 2|4) in different ways and must be studied separately. 

As in the two-impurity case, the problem of diagonalizing the perturbation simplifies 
enormously when the matrix elements are expanded to leading order in A'. We will take 
this approach here to obtain an overview of how degeneracies are lifted by the interaction. 
The generalization of the results to all loop orders in A' (but still to first non-leading order 
in 1/J) will be presented in the next section. We use the same methods, conventions and 
notations as in our recent detailed study of the two-impurity problem |2] (especially in Sec. 6 
of that paper). It is once again the case that in the one- loop approximation, projection onto 
invariant subspaces under the manifest global SO (A) x SO (4) symmetry often diagonalizes 
the Hamiltonian directly (and at worst reduces it to a low-dimensional matrix). Symbolic 
manipulation programs were used to organize the complicated algebra and to perform explicit 
projections onto invariant subspaces. 
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3.1 Matrix evaluation: inequivalent mode indices (q ^ r ^ s) 

In the sector of spacetime bosons, the subspace of purely bosonic states a^a^a^ \J) is 
5 12- dimensional. When each of the three mode indices (q,r, s) are different, states with bi- 
fermionic excitations a^b^b^ |J) are inequivalent under permutation of the mode indices, 
and form a 1,536-dimensional subsector. The entire bosonic sector of the three-impurity 
state space therefore contains 2,048 linearly independent states. The fermionic sector de- 
composes in a similar manner: the subsector of purely fermionic states b^bybj* | J) is 512- 
dimensional; fermionic states containing two bosonic excitations a^a^b^ \J) are inequiva- 
lent under permutation of the mode indices, and comprise an additional 1,536-dimensional 
subsector. Adding this 2,048-dimensional fermion sector brings the dimensionality of the 
entire state space to 4,096. 

Our first task is to evaluate the interaction Hamiltonian matrix. The matrix elements 
needed to fill out the spacetime boson sector are listed in table 01 To evaluate the entries, we 
express the Hamiltonians (j2.15H2.l7jl in terms of mode creation and annihilation operators, 
expand the result in powers of A' and then compute the indicated matrix elements between 
three- impurity Fock space states. We collect below all the relevant results of this exercise. 





af a faft|j) 


a?b?b$\J) 


a^b s J\J) 


a^bfb^J) 


(J\ afafaC 


H-qb 


Hbf 


Hbf 


Hbf 




Hbf 


Hff + Hbf 


Hbf 


Hbf 


(J\afbffl 


Hbf 


Hbf 


Hff + Hbf 


Hbf 


{J\aXK 


Hbf 


Hbf 


Hbf 


Hff + Hbf 



Table 3: Interaction Hamiltonian on spacetime boson three-impurity string states (g / r ^ s) 



We will use an obvious (m, n) matrix notation to distinguish the different entries in 
table El The purely bosonic, 512-dimensional (1, 1) block has the explicit form 



(j\aX^ (H BB )a?a?a?\ J) = ^-5 AF 5 BE 5 CD (rs + q(r + a) - g 2 

'J 



2 ) 



+ — l) Ar 
2J 



+2rs(5 bd 5 ce - 6 bc 5 de -5 b ' d '5 c ' e ' + 5 b ' c '5 d ' e ')] + (r ^ q, F ^ E, A ^ B) 

+ ( a F^D, A^C)\ . (3.3) 



The off-diagonal entries that mix purely bosonic states j' | J) with states containing 

bi-fermions a^b^b^ \J) are given by a separate set of 512-dimensional matrices. The (1, 2) 
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block in table EJ for example, yields 



a'b 



7 



a/3 



(3.4) 



where the index a (a') symbolizes the value of the vector index A, provided it is in the first 
(second) SO (A). There are six blocks in this subsector, each given by a simple permutation 
of the mode indices (g, r, s) in eqn. ()3.4j) . In table El these matrices occupy the (1, 2), (1, 3) 
and (1,4) blocks, along with their transposes in the (2, 1), (3, 1) and (4, 1) entries. 

The pure-fermion sector of the Hamiltonian, H-p-p, has non- vanishing matrix elements 
between states containing bi-fermionic excitations. The H FF contribution to the (2, 2) block, 
for example, is given by 



[j\bffiaf(H FF )a^b?b 5 J\j) 



2J 



(q - r) 2 6 AB 6 aS 6^ 



fr*) " - - (7TV) 



24J 



Pi 



a/3 . 



(3.5) 



A similar contribution, related to this one by simple permutations of the mode indices (q, r, s), 
appears in the diagonal blocks (3,3) and (4,4) as well. 

The bose-fermi mixing Hamiltonian H BF makes the following contribution to the lower 
diagonal blocks (2,2), (3,3) and (4,4) in table EJ 



A' 



[j\bffiaf(H BF )af%?b s J\j) = ^\2s(q + r-s)5 ab 5 a5 5^ - rs ( 7 ab ) Pl - ( 7 a ' fe ') 



-sq 



(7°T - ( r b ) 



q 2 + r 2 + s 2 — s(q + r) 



S a ' b 'S a5 S^ 



(3.6) 



The i?BF sector also makes the following contribution to the off-diagonal (2, 3) block: 

(j\bffla A (H BF )a B ^b s J \j) = -^° 7 "{ - V*")** - {l ab f 5 + (7 a ' b T} ■ (3-7) 

The contributions of Hbf to the remaining off-diagonal blocks (2, 3), (2, 4), etc. are obtained 
by appropriate index permutations. 

The sector of spacetime fermions decomposes in a similar fashion. The fermion analogue 
of tableElfor the bosonic sector appears in tableEJ The (1, 1) fermion block is occupied by the 
pure-fermion sector of the Hamiltonian taken between the purely fermionic three-impurity 
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Hint 


b?b*b$\j) 






°?*a?l#\J) 


(J\ b°b?b] 


H FF 




Hbf 


Hbf 


( J\ ftjafrf 


Hbf 


Hbb + Hbf 


Hbf 


Hbf 


(J|6>X 


Hbf 


Hbf 


Hbb + Hbf 


Hbf 


(J\b?a A a B 


Hbf 


Hbf 


Hbf 


Hbb + Hbf 



Table 4: Interaction Hamiltonian on spacetime fermion three-impurity states (g ^ r / s) 
states bfb^bf \J): 



(j\b a q b?b:(H FF )bi^b 5 J\j) 
A' 



^ [q 2 + r 2 + s 2 -rs- q{r + s)] S a5 8^ 



+ ( r ^ g, a?i|3, 5 ^ e) + (s ^ g, a ^ 7, 5 ^ C)| • (3-£ 



The off-diagonal (1,2), (1,3) and (1,4) blocks (and their transposes) mix purely fermionic 
states with a A ^a B ^b°^ \ J) states: 



[j\b«b?b:(H BF )afa^\j) 



~ 7 



,a'b 



01 



(3.9) 



The lower-diagonal (2,2), (3,3) and (4,4) blocks receive contributions from the pure boson 
sector of the Hamiltonian: 

j\b°a?a?(H BB )a?a?bP\j) = -^j^Ut - s) 2 5 BC 5 AD - (r 2 + s 2 )(d ad 5 bc - <W C ') 



-2rs (d ac 5 bd - 5 ab 5 cd -5 a ' c '5 b ' d ' + 5 a ' b '5 c ' d '^ + (r 2 - s 2 ) (d ad 5 

In the same diagonal blocks of table 0J the Hbf sector contributes 

A' 



2^ / xadxb'c' ^a'd' ^bc 



.(3.10) 



zBC 



(j\b«a?af(H BF )a?a?^\j) = [+(Bq(r + s) - 5(r 2 + s 2 ) - 6q 2 )5 AD 5 

+ (3q 2 + s 2 )6 AD 5 bc +(3q 2 + r 2 )6 BC 5 ad + (r 2 - 5q 2 )6 BC 5 a ' d ' + (s 2 - hq 2 )5 AD 5 V ' 

-A5 BC qr\ (r d ) a(3 - {r' d T P } -4^g S [( 7 6c ) Q/3 - ( 7 bV ) a/3 \ . (3.11) 
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Finally, the off-diagonal blocks (2, 3), (2, 4) and (3, 4) (plus their transpose entries) are given 
by the H B -p matrix element 



[j\b a q a A af(H BF )a^af^\j) 



X 



5 AC\ 5 a/3 



32J 

-(q 2 -I8qs + s 2 )5 b ' d ' 



+ lQqs 



s f 6 BD _ (g 2 + Ugs + S 2 )5 M 
bd\ a P („,b'd'\ a P 



U>dyp _ ^b'd'y 



(3.12) 



A significant departure from the two- impurity case is that all these matrix elements have, 
along with their spacetime index structures, non-trivial dependence on the mode indices. The 
eigenvalues could potentially have very complicated mode-index dependence but, as we shall 
see, they do not. This amounts to a rigid consistency check on the whole procedure that 
was not present in the two-impurity case. 



3.2 Matrix diagonalization: inequivalent mode indices (q ^ r ^ s) 

We now turn to the task of diagonalizing the one-loop approximation to the perturbing 
Hamiltonian. To simplify the task, we exploit certain block diagonalizations that hold to 
leading order in A' (but not to higher orders). While we eventually want to study the spec- 
trum to all orders in A', diagonalizing the Hamiltonian at one loop will reveal the underlying 
supermultiplet structure. As an example of the simplifications we have in mind, we infer 
from ()3.4|) that the matrix elements of H BF between pure boson states a^a^a^ \J) and 
bifermionic spacetime bosons vanish to leading order in A' if all three 5*0(8) bosonic vector 
indices lie within the same SO (4), descended either from AdS^ or S 5 . Restricting to such 
states brings the bosonic sector of the Hamiltonian into the block-diagonal form in table El 
This leaves two 64-dimensional subspaces of purely bosonic states on which the perturbation 





a^aPa* | J) + a a 'ta b V't | J) a A ^b a ^ \ J) 


(J\a a a b a c + (J\a a 'a b 'a c ' 
(J\ a A b a W 


H BB 
Hff + H B f 



Table 5: Block-diagonal S0(4) projection on bosonic three-impurity string states 
is block diagonal, as recorded in table [HI 



Hint 


a at a b t a c t | J) a'W 1 " \J) 


(J\ a a a b a c 
(J\ a a 'a b 'a c ' 


(^Bb)64x64 

(-f^BB)64x64 



Table 6: SO (4) projection on purely bosonic states 
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Since the interaction Hamiltonian has manifest SO (4) x SO (4) symmetry, it is useful to 
project matrix elements onto irreps of that group before diagonalizing. In some cases the irrep 
is unique, and projection directly identifies the corresponding eigenvalue. In the cases where 
an irrep has multiple occurrences, there emerges an unavoidable matrix diagonalization that 
is typically of low dimension. In what follows, we will collect the results of carrying out this 
program on the one-loop interaction Hamiltonian. A very important feature of the results 
which appear is that all the eigenvalues turn out to have a common simple dependence on 
mode indices. More precisely, the expansion of the eigenvalues for inequivalent mode indices 
(q, r, s) out to first non-leading order in A' and 1/J can be written as 

r, .) = 3 + ^V±ii±f!) ^ + a + Q( j. 2) j (3 13) 

where A is a pure number that characterizes the lifting of the degeneracy in the various 
sectors. The notation Abb, Abf and Aff will be used to denote energy corrections arising 
entirely from the indicated sectors of the perturbing Hamiltonian. This simple quadratic 
dependence of the eigenvalues on the mode indices does not automatically follow from the 
structure of the matrix elements themselves, but is important for the successful match to 
gauge theory eigenvalues. In what follows, we will catalog some of the different A values 
that occur, along with their SO (4) x 5*0(4) irreps (and multiplicities). When we have the 
complete list, we will discuss how they are organized into supermultiplets. 

In the 5*0(4) projection in table El we will find a set of 64 eigenvalues for both the 
50(4)^ d 5 and 50(4) 5 s subsectors. We record this eigenvalue spectrum in table using 
an SU{2) 2 x SU{2) 2 notation. For comparison, it is displayed alongside the projection 
of the 2-impurity spectrum onto the same subspace (as found in [HJ E]). In the three- 



SO(A) AdS x 50(4) 55 


Abb 




SO(A) AdS x 50(4) 55 


Abb 


(1,1; 2, 2) 


-8 




(i,i; i,i) 


-6 


[1,1; (2 + 4), 2] + [1,1; 2, (2 + 4)] 


-6 




(1,1; 3,1) + (1,1; 1,3) 


-4 


[1,1; (2 + 4), (2 + 4)] 


-4 




(i,i; 3, 3) 


-2 


[(2 + 4), (2 + 4); 1,1] 


-2 




(3, 3; 1,1) 


-2 


[(2 + 4), 2; 1,1] + [2, (2 + 4); 1,1] 







(3,1; 1,1) + (1,3; 1,1) 





(2, 2; 1,1) 


2 




(i,i; i,i) 


2 



Table 7: Three-impurity energy spectrum in the pure-boson 50(4) projection (left panel) 
and two-impurity energy spectrum in the same projection (right panel) 

impurity case, the (1, 1; 2, 2) level in the 50(4) 5 s subsector clearly descends from the two- 
impurity singlet (1, 1; 1, 1) in the same 50(4) subgroup. In the same manner, the three- 
impurity [1, 1; (2 + 4), 2] + [1, 1; 2, (2 + 4)] level descends from the 50(4)^5 antisymmetric 
two-impurity state (1, 1; 3, 1) + (1, 1; 1, 3), and the three- impurity [1, 1; (2 + 4), (2 + 4)] 
level is tied to the two-impurity symmetric-traceless (1,1; 3, 3) irrep. In the 50(4) 5 5 sub- 
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sector, each of these levels receives a shift to the energy of —2. The total multiplicity of 
each of these levels is also increased by a factor of four when the additional (2, 2) is ten- 
sored into the two- impurity state space. The SO(4)Ads subsector follows a similar pattern: 
the (2, 2; 1, 1), [(2 + 4), 2; 1, 1] + [2, (2 + 4); 1, 1] and [(2 + 4), (2 + 4); 1, 1] levels appear 
as three-impurity descendants of the two-impurity irrep spectrum (1, 1; 1, 1) + (3, 1; 1, 1) + 
(1, 3; 1, 1) + (3, 3; 1, 1). In this subsector, however, the three-impurity energies are identical 
to those in the two-impurity theory. 

The bosonic 50(4) projection has a precise fermionic analogue. Similar to the bosons, 
the £0(9, 1) spinors b q decompose as (2, 1; 2, 1) + (1, 2; 1, 2) under the action of II parity: 

n&J = i\ m\ = -l\ . (3.i4) 

The notation b q labels (1, 2; 1, 2) spinors with positive eigenvalue under II, and b q indicate 
(2, 1; 2, 1) spinors which are negative under II. Analogous to the £0(4) projection on the 
£0(8) bosonic operators a A ^ — > + a q t, projecting out the positive or negative eigenvalues 
of II on the eight-component spinor leaves a subspace of four- component spinors spanned 

by b q and b q . 

We can perform a projection on the subsector in table El similar to that appearing in 
table El In this case, instead of three bosonic impurities mixing with a single bosonic (plus 
a bi-fermionic) excitation, we are now interested in projecting out particular interactions 
between a purely fermionic state and a state with one fermionic and two bosonic excitations. 
Using ± to denote the particular representation of the fermionic excitations, the off-diagonal 

elements given by ()3.9|) vanish for + + H > ± and > ± interactions. In other words, 

the pure fermion states in the (1,1) block of table H] will not mix with states containing 
two bosonic excitations if all three fermionic oscillators lie in the same II projection. This 
projection appears schematically in table El 



Hmt 


y*\yt\yt\ \ j) + \ j) a A ^a B ^ \ j) 


(J\ b a Pfr + (J\ b a ¥w 
(J\ a A a B b a 


H FF 
H BB + H BF 



Table 8: Block-diagonal projection on fermionic three-impurity string states 

The (1,1) pure fermion block in table El breaks into two 64-dimensional subsectors under 
this projection. By tensoring an additional (1, 2; 1, 2) or (2, 1; 2, 1) impurity into the two- 
impurity state space, we expect to see a multiplicity structure in this projection given by 

(1,2) x (1,2; 1,2) = (1,2; 1,2) + [1,2; 1,(2 + 4)] 

+ [1,(2 + 4); 1,2] + [1,(2 + 4); 1,(2 + 4)] , 

(2,1) x (2,1; 2,1) = (2,1;2,1) + [2,1;(2 + 4),1] 

+ [(2 + 4), 1; 2, 1] + [(2 + 4), 1; (2 + 4), 1] , (3.15) 
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for a total of 128 states. The projections onto the two 64-dimensional n + and Il_ subspaces 
yield identical eigenvalues and multiplicities. The results for both subspaces are presented 
in table EJ The two-impurity bi-fermion states in table El are spacetime bosons while the 



SO(4) AdS x SO(4) 55 


Aff 




(2,1; 2,1) + (1,2; 1,2) 


-3 


[2,1; 


(2 + 4), 1] + [1,2; 1,(2 + 4)] 


-1 


[(2 + 


4), 1; 2,1] + [1,(2 + 4); 1,2] 


-5 


[(2 + 4),l; 


(2 + 4), 1] + [1,(2 + 4); 1,(2 + 4)] 


-3 



SO(A) AdS x 


SO(4) s5 


A FF 


(1,1; 1,1) + 


(i,i; 1,1) 


-2 


(1,1; 3,1) + 


(1,1; 1,3) 





(3,1; 1,1) + 


(1,3; 1,1) 


-4 


(3,1; 3,1) + 


(1,3; 1,3) 


-2 



Table 9: Spectrum of three-impurity states (left panel) and two-impurity states (right panel) 
created by n_|--projected fermionic creation operators 

tri-fermion states are spacetime fermions. For comparison purposes, we have displayed both 
spectra. Note that the 0(1/ J) energy corrections of the two types of state are simply 
displaced by —1 relative to each other. 

This exhausts the subspaces that can be diagonalized by simple irrep projections. The 
remaining eigenvalues must be obtained by explicit diagonalization of finite dimensional 
submatrices obtained by projection onto representations with multiple occurrence. The 
upshot of these more complicated eigenvalue calculations is that the first-order A' eigenvalues 
take on all integer values from A = —8 to A = +2, alternating between spacetime bosons 
and fermions as A is successively incremented by one unit. 



3.3 Assembling eigenvalues into supermultiplets 

Finally, we need to understand how the perturbed three-impurity spectrum breaks up into 
extended supersymmetry multiplets. This is relatively easy to infer from the multiplicities 
of the perturbed eigenvalues (and the multiplicities are a side result of the calculation of the 
eigenvalues themselves). In the last subsection, we described a procedure for diagonalizing 
the one-loop perturbing Hamiltonian on the 4, 096-dimensional space of three-impurity string 
states with mode indices p ^ q ^ r. The complete results for the eigenvalues A and their 
multiplicities are stated in table El (we use the notation of ()3.13|) . while the B and F 
subscripts are used to indicate bosonic and fermionic levels in the supermultiplet). 



A 


-8 -7 -6 -5 -4 -3 -2 -10 12 


Multiplicity 


4 B 40 F 180 B 480 F 840 B 1008 F 840 B 480 F 180 B 40 F 4 B 



Table 10: Complete three- impurity energy spectrum (with multiplicities) 
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The eigenvalues A must meet certain conditions if the requirements of PSU(2, 2|4) sym- 
metry are to be met. The eigenvalues in question are lightcone energies and thus dual to 
the gauge theory quantity A = D — J, the difference between scaling dimension and R- 
charge. Since conformal invariance is part of the full symmetry group, states are organized 
into conformal multiplets built on conformal primaries. A supermultiplet will contain several 
conformal primaries having the same value of A and transforming into each other under the 
supercharges. All 16 supercharges increment the dimension of an operator by 1/2, but only 
8 of them (call them Q a ) also increment the .R-charge by 1/2, so as to leave A unchanged. 
These 8 supercharges act as 'raising operators' on the conformal primaries of a supermulti- 
plet: starting from a super-primary of lowest -R-charge, the other conformal primaries are 
created by acting on it in all possible ways with the eight Q a . Primaries obtained by acting 
with L factors of Q a on the super-primary are said to be at level L in the supermultiplet 
(since the Q a anticommute, the range is L = to L = 8). The multiplicities of states at 
the various levels are also determined: for each L = primary operator, there will be Of 
such operators at level L (where O^ is the binomial coefficient). If the L = primary has 
multiplicity s, summing over all L gives 2 8 s = 256s conformal primaries in all. 

These facts severely restrict the quantity A in the general expression ()3.13p above. Al- 
though the states in the degenerate multiplet all have the same J, they actually belong to 
different levels L in more than one supermultiplet. A state of given L is a member of a 
supermultiplet built on a 'highest- weight' or super-primary state with R = J — L/2. Since 
all the primaries in a supermultiplet have the same A, the joint dependence of eigenvalues 
on A, J,L must be of the form A(A, J — L/2). The only way the expansion of (|3.13|) can 
be consistent with this is if A = L + c, where c is a pure numerical constant (recall that 
A' = A/ J 2 ). Successive members of a supermultiplet must therefore have eigenvalues sepa- 
rated by exactly 1 and the difference between 'top' (L = 8) and 'bottom' (L = 0) eigenvalues 
for A must be exactly 8. 

The A eigenvalues in table El are integer-spaced, which is consistent with supersymmetry 
requirements. However, because the range between top and bottom eigenvalues is 10, rather 
than 8, the 4, 096-dimensional space must be built on more than one type of extended 
supermultiplet, with more than one choice of c in the general formula A = L + c. This is 
to be contrasted with the two-impurity case, where the degenerate space was exactly 256- 
dimensional and was spanned by a single super conformal primary whose lowest member was 
a singlet under both Lorentz transformations and the residual SO (4) i?-symmetry. We can 
readily infer what super conformal primaries are needed to span the degenerate three-impurity 
state space by applying a little numerology to table EH The lowest eigenvalue is A = —8: it 
has multiplicity 4 and, according to table its 5*0(4) x SO (A) decomposition is (1, 1; 2, 2) 
(spacetime scalar, .R-charge 5*0(4) fo ur- vector) . According to the general arguments about 
how the full extended supermultiplet is built by acting on a 'bottom' state with the eight 
raising operators, it is the base of a supermultiplet of 4 x 256 states extending up to A = 0. By 
the same token, there is a highest eigenvalue A = +2: it has multiplicity 4 and, according to 
table0 its 50(4) x 50(4) decomposition is (2, 2; 1, 1) (spacetime vector, Recharge singlet). 
Using lowering operators instead of raising operators, we see that one derives from it a 
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supermultiplet of 4 x 256 operators with eigenvalues extending from A = —6 to A = +2. 
The multiplicities of the A eigenvalues occurring in these two supermultiplets are of course 
given by binomial coefficients, as described above. By comparing with the total multiplicities 
of each allowed A (as listed in table ITUj) we readily see that what remains are 8 x 256 states 
with eigenvalues running from A = —7 to A = +1 with the correct binomial coefficient 
pattern of multiplicities. The top and bottom states here are spacetime fermions and must 
lie in a spinor representation of the Lorentz group. It is not hard to see that they lie in the 
eight-dimensional 50(4) x SO (A) irrep (2, 1; 1, 2) + (1, 2; 2, 1). This exhausts all the states 
and we conclude that the three-impurity state space is spanned by three distinct extended 
superconformal multiplets. The detailed spectrum is given in table ^2 (where the last line 
records the total multiplicity at each level as given in table [HI and the first line records 
the two-impurity spectrum for reference). Note the peculiar feature that certain energies 
are shared by all three multiplets: this is an accidental degeneracy that does not survive at 
higher loop order. 



A 


-8 


-7 


-6 


-5 


-4 


-3 


-2 


-1 





1 


2 




A = 2 






Ib 


8f 


28 B 


56 F 


70 B 


56 F 


28 B 


8 F 


Ib 


scalar 


A = 3 


4 


32 


112 


224 


280 


224 


112 


32 


4 






6*0(4)55 vector 








4 


32 


112 


224 


280 


224 


112 


32 


4 


SO(4:) AdS5 vector 






8 


64 


224 


448 


560 


448 


224 


64 


8 




spinor 


Total 


4 B 


40 F 


180 B 


480 F 


840 B 


1008 F 


840 B 


480 F 


180 B 


40 F 


4b 


4,096 



Table 11: Submultiplet breakup of the three- impurity spectrum 

A complete analysis of the agreement with gauge theory anomalous dimensions will have 
to be deferred until a later section: the dimensions of three-impurity gauge theory operators 
are much harder to calculate than those of the two-impurity operators and there are few 
results in the literature, even at one loop. However, it is worth making a few preliminary 
remarks at this point. Since there are three superconformal multiplets, we have only three 
independent anomalous dimensions to compute. Minahan and Zarembo [3] found that the 
problem simplifies dramatically if we study the one-loop anomalous dimension of the special 
subset of single-trace operators of the form tr ((f) 1 Z J ) (and all possible permutations of the 
fields inside the trace), where the -R-charge is carried by an 50(4) x 50(4) singlet scalar 
field Z and the impurities are insertions of a scalar field (f) lying in the (1, 1; 2, 2) (vector) 
irrep of the residual 50(4) x 50(4) symmetry. More formally, these operators are in the 
50(4) x 50(4) irrep obtained by completely symmetrizing / vectors in the (1,1; 2, 2) ir- 
rep. The crucial point is that such operators form a 'closed sector', mixing only among 
themselves under the anomalous dimension operator. More importantly, the action of the 
one-loop anomalous dimension operator on this closed sector can be recast as the action of an 
integrable spin chain Hamiltonian of a type solvable by Bethe ansatz techniques. Although 
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the Bethe ansatz is generally not analytically soluble, Minahan and Zarembo used it to ob- 
tain a virial expansion for the anomalous dimension in which the number I of impurities is 
held fixed, while the i?-charge J is taken to be large (see eqn. (5.29) in In terms of the 
number of spin chain lattice sites £, their result appears as 

7*o(6) = ^ Y, Mnk " (£ + M n + l) + 0(£" 4 ) . (3.16) 

n 

The integer k n represents pseudoparticle momenta on the spin chain, and is dual to the string 
theory worldsheet mode indices; the quantity M n labels the number of trace impurities with 
identical k n . With / impurities, the spin chain length is given in terms of the i?-charge by 
£ = J + I, which leads to 

7*(6) = E Mnk " ( J ~ 2I + M n + 1) + 0(J~ 4 ) . (3.17) 

n 

This virial expansion is similar in character to 1)3.13)1 and, for 1 = 3 (the three-impurity 
case), it matches that equation precisely with A = —4. 

On the string theory side, three completely symmetrized (1, 1; 2, 2) vectors form a tensor 
in the (1, 1; 4, 4) irrep; such an irrep can be constructed from three SO (4)^5 vector (bosonic) 
creation operators. Table [7| shows that the corresponding string perturbation theory eigen- 
value is (at one-loop order) A = —4 as well. We infer from table ^2 that this eigenvalue lies 
at level L = 4 of the 5*0(4)55 vector superconformal multiplet (and this argument takes care 
of the gauge theory/string theory comparison for all other operators in that multiplet). 

The sector described above is often called an so(6) 2 sector on the gauge theory side, with 
reference to the subalgebra of the full superconformal algebra under which it is invariant. 
In an su(2) subspace of the so (6), this sector becomes closed to all loop order. For future 
reference, we note that Beisert [12^ has identified two other 'closed sectors' of operators in the 
gauge theory. In addition to the bosonic su(2) sector, a bosonic $1(2) sector and an su(2|3) 
sector (of which the closed su(2) sector is a subgroup) are also exactly closed. It should be 
noted that integrable $1(2) spin chains were discovered some time ago in phenomenologically- 
motivated studies of the scaling behavior of high-energy scattering amplitudes in physical, 
non-supersymmetric QCD [T3] (see also [T^l [TBI [XZ| ) - The su(2|3) spin chain was studied 
more recently in 6j: this closed sector breaks into the $u(2) bosonic sector and a fermionic 
subsector which, to avoid confusion, we simply denote as a subgroup of su(2|3). 

In the string theory, the subsectors analogous to the gauge theory $1(2) and su(2|3) are 
constructed out of completely symmetrized SO(4)acls bosons and completely symmetrized 
fermions of the same II eigenvalue, respectively They correspond to the central L = 4 
levels of the remaining two supermultiplets in table ^2 an d a calculation of their eigenvalues 
would complete the analysis of the match between three-impurity operators and string states 

2 This notation is used to distinguish the protected gauge theory symmetry groups from those in the string 
theory. 
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at one-loop order. Unfortunately, explicit general results for three-impurity operator dimen- 
sions, analogous to those obtained by Minahan and Zarembo for the 50(6) sector [3], have 
not been obtained for the other closed sectors. The Bethe ansatz for the general one-loop 
integrable spin chain presented in jXBj could easily be exploited for this purpose. 3 However, 
since we eventually want to go beyond one-loop, where Bethe ansatz technology is less well- 
developed, we have found it more useful to develop numerical methods for evaluating spin 
chain eigenvalues (we refer the reader to [IlJ for a check of our results against Bethe- ansatz 
techniques, including the higher-loop corrections of ^H])- This subject will be developed in 
a later section. 

3.4 Two equivalent mode indices (q = r = n, s = — 2n) 

When two mode indices are allowed to be equal, the analysis becomes slightly more com- 
plicated. Since we are diagonalizing a Hamiltonian that is quartic in oscillators in a basis 
of three-impurity string states, one oscillator in the "in" state must always be directly con- 
tracted with one oscillator in the "out" state and, with two equal mode indices, there are 
many more nonvanishing contributions to each matrix element. While the matrix elements 
are more complicated, the state space is only half as large when two mode indices are allowed 
to be equal (only half as many mode-index permutations on the basis states generate lin- 
early independent states). As a result, the fermionic and bosonic sectors of the Hamiltonian 
are are each 1,024-dimensional. By the same token, the multiplet structure of the energy 
eigenstates will be significantly different from the unequal mode index case studied in the 
previous subsection. 

To study this case, we make the mode index choice 

q = r = n s = —In . (3.18) 

The structure of matrix elements of the string Hamiltonian between spacetime bosons is 
given in table ED This table seems to describe a 3 x 3 block matrix with 512 x 512 blocks 



#hit 


a-\ n a E Ja F J\J) a D _ljPX\J) aXb% n \J) 


(J a^a I ^a c _ 2n 


Hbb Hbf Hbf 
Hbf Hff + Hbf Hbf 
Hbf Hbf Hff + Hbf 



Table 12: Bosonic three-impurity string perturbation matrix with (q = r = n, s = —2n) 

in each subsector, giving a 1,536-dimensional state space. However, the vector and spinor 
indices are required to run over values that generate linearly independent basis states. This 

3 We are indebted to N. Beisert for this observation. 
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eliminates one third of the possible index assignments, implying that the matrix is in fact 
1,024 x 1,024. 

To evaluate the entries in table we express the Hamiltonians (|2. 15112. lTJ) in terms of 
mode creation and annihilation operators, expand the result in powers of A' and compute 
the indicated matrix elements between three-impurity Fock space states. We collect below 
all the relevant results of this exercise for this equal-mode-index case. 

The purely bosonic subsector in the (1,1) block is given by 



(J\a A a B a% n (H BB )a^l n a^a F1 \J) 



n 2 X 
~2l 



5 S BF S cd S ae + 5 S AF S cd 5 be - 4 b Bt S ad S 



AF red i 



:BF zadzce 



+4 6 BF 5 ac 5 de + 4 5 AF 5 bc 5 de + 5 5 BE 5 cd 5 af - 4 5 BE 5 ad 5 cf + 4 5 BE 5 ac 5 df + 4 6 AE 5 bc t 
-4 S bd (V F S ce + 5 AE 5 C/ ) + 3 S BF 5 ae S c ' d ' + 3 S AF S be S c ' d ' + 3 S BE 5 af S c ' d ' 

-3 5 BF 5 cd 5 a ' e ' - 3 5 AF 5 cd 5 b ' e ' - 5 5 BF 5 c ' d 'd a ' e ' - 5 5 AF 6 c ' d 'd b ' e ' + 4 6 BF 5 a ' d '5 c ' e ' 
+4 S AF 5 b ' d 'S c ' e ' - 4 5 BF S a ' c 'S d ' e ' - 4 S AF S b ' c 'S d ' e ' - 3 S BE S cd 5 a ' f ' - 3 S AE S cd S b ' f ' 
-5 6 BE 5 c ' d '5 a 'f - 5 5 AE 5 c ' d '5 b 'f + 4 5 BE 5 a ' d '5^' + 4 6 AE 6 b ' d '5^' - 4 6 BE 5 a ' c '5 d 'f 
-4 S AE S b ' c 'S d ' f ' + S AE S bf (5 S cd + 3 5 c ' d ') - 2 S CD h ( S BE S AF + S AE 5 BF ) - 5 be S af 

-5 ae 6 bf + 5 ab 5 ef + 5 b ' e '5 a ' f ' + 5 a ' e '5 b ' f ' - 5 a ' b '5 e ' f 



(3.19) 



This matrix element exhibits the same antisymmetry between the SO(4)Ads and SO(4)s5 
indices that is exhibited in eqn. ()3.3|) . The off-diagonal H BF mixing sector is essentially 
equivalent to its counterpart in eqn. (|3.4|) : 



J 



a n a n ci_2 n 



(H BF )a%bX 



J 



n 2 \' f / A a/3 



2J 



7 



a'b 



7 



a/3 



(3.20) 



The diagonal contributions from the pure fermion sector Hpp m the (2, 2) and (3, 3) blocks 
of table H21 appear as 



J 



Kb^ n (H FF )a B _lbX 



J 



n 2 X' 



24J 



S AB ) (fly 7 (ry^y 5 - ( 7 «) a/3 ( 7 «) 7 ' 5 



_( 7 ii) Q<5 (yi)^_ (//)^(//)^V(/iy^(//)^ + (/iy o (/iy^j .(3.21) 

The H BF sector exhibits the following contribution to the lower diagonal blocks (2, 2) and 
(3,3): 



a/3 



\ 7<* 



J 



««- 2 ,(#BF)afL« 



J 



n 2 X ' 



-10 5 a ' b ' (5 a5 5^ - 5 a ~<5 p& ) 



-8 5 ab (5 a5 5^ - 5 ai 5 p& ) -5 ai (j ab ) f3S - ( 7 a ' 6 ') 



\/3S 



+ 5' 



CtS 



+6* 



( 7 ab ^7 _ ( 7 a ' b ') a 7 



(3.22) 
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Finally, the off-diagonal version of ()3.22|) appears in the (2, 3) block (along with its transpose 
in the (3,2) block): 

n 2 \' 



J 



b a n b^ 2n (H BF )aXb% n 



J 



J 



fia'b' (fiOtS^Pi _ fi a lfiP s ^ _ fiab f^aS ^7 _ ja7£/3<5\ 



-5^ 



^abyS _ ^a'b'yS 



(3.23) 



The fermionic sector perturbation matrix is displayed schematically in table ED Like 
table ^1 it is 1,024 x 1,024 once redundant index assignments are eliminated. 



Hint 


b% n bX\J) a c \a^\J) a c Ja^b% n \J) 


(J\b a Xbl 2n 
(J\ b a _ 2 M 


Hff Hbf Hbf 
Hbf Hbb + Hbf Hbf 
Hbf Hbf Hbb + Hbf 



Table 13: Fermionic three-impurity string perturbation matrix 



r = n, s = —2n) 



The purely fermionic subsector in the (1, 1) block of table El takes the form 

9n 2 X 



J 



K€bl 2n {HFF)b% n bX 



J 



J 



-<5 7C (5 ae 5 ps - 8 a5 8 pe ) 



n 2 \' 
24J 



(ryV\<*P(ryijyt _ fyti ) ocS /3e _|_ /^ijye/yijyS _ f^i'j'yP^i'j'yt 



25 



a.5 



+ fyijyC(ryijye _(yi'j'yyfyi'fyC _|_ ^yi'j'y^yi'j'yC _ (7*' J '')^(7*'-'') 



+25 c 



(jvy-Y^yijyc -(jijysfyiiys + (yvyCfyvys _ (yvryi^ytr^ 



28> 3S 



fyijyyfyijyC _ fyijyef^ijyt 



+ fyijyCfyijye _fyi'j'\orf/yi'j'yC + fyi'j'yef^i'j'yC _ fyi'j'yCfyi'j'ye 



-28* 



fyijyjfyijyc _/yijyStyijyC + /yijy(/yijyS _ /yi'j'\ory/yt'3'yc 



The off-diagonal blocks (1,2) and (1,3) receive contributions from the Hbf sector: 



(3.24) 



J 



b a Xbl 2n {HBF)a A _\ n a B n ^ 



J 



J 



a 6 



aJb' 



7 



y ab ' 

cry 



a'b 



7 



y a ' b 
ay 



Pr 



(3.25) 
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The bosonic sector H BB contributes to the (2, 2) and (3, 3) blocks: 

n 2 \' 



cd 



{j\b a q 4a*{H BB )afa°X\j) = --^j-5 a ^9 6 AD 6 BC + A5 ac 6 bd - A5 ab 6' 

_ 6 ad ^ 6 bc + 3 6 b'c'^ _ 4 6 a'c' 6 b'd' + 4 S a'b' 6 c'd> + ga'd' ^ £&V + 3 gbc^ 1 _ (33^ 

In the same lower-diagonal blocks, H B f exhibits the contribution 

n 2 X' 



J 



by*a* 2n (H BF )a c _laX 



J 



8J 



S. 39 5^5 AD 5 BC + rt AD (> v - 75 bc 



_ A5 aP S BC 6 ad _ 5 a>d> , +4 tf 



-85 



AD 



(3.27) 



Finally, H B f yields matrix elements in the off-diagonal block (2,3): 

n 2 \' 



J 



by A a B 2n (H BF )a^a^_\ n 



J 



32J 



QS aP 5 AC S BD + Q§ af3 S AD s BC 



+6 a/3 6 AC ^bd _ A16 b'd'^ + 5 aP 5 AD ^bc _ ^b'c' 



-325 



AD 



325 



AC 



(3.28) 



We can perform a full symbolic diagonalization of the 1, 024 x 1, 024 bosonic and fermionic 
perturbation matrices to obtain the one-loop in A', 0(1/ J) energy corrections. They can 
all be expressed in terms of dimensionless eigenvalues A according to the standard formula 
()3.13|) modified by setting q = r = n, s = — 2n: 



E,(n) = • 3n 2 A' [ \ ■ 1 .()(./ 2 



(3.29) 



The resulting spectrum is displayed in table The levels clearly organize themselves into 



Ai (S 5 vector) 


-23/3 


-20/3 


-17/3 


-14/3 


-11/3 


-8/3 


-5/3 


-2/3 


1/3 


Multiplicity 


4b 


32 F 


112 B 


224 F 


280 B 


224 F 


112 B 


32b 


4b 


A 2 (AdSs vector) 


-19/3 


-16/3 


-13/3 


-10/3 


-7/3 


-4/3 


-1/3 


2/3 


5/3 


Multiplicity 


4b 


32 F 


112b 


224 F 


280 B 


224 F 


112b 


32b 


4b 



Table 14: Spectrum of three-impurity string Hamiltonian with (q = r = n, s = —2n) 
two superconformal multiplets built on vector primary states. Note that the spinor multiplet 
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is absent and that the degeneracy between multiplets that was seen in the inequivalent 
mode index case has been lifted. The spinor multiplet is absent for the following reason: it 
contains a representation at level L = 4 arising from fermion creation operators completely 
symmetrized on 5*0(4) x 5*0(4) spinor indices; such a construct must vanish unless all the 
creation operator mode indices are different. 

If we keep track of the 50(4) x 50(4) irrep structure, we find that the symmetric-traceless 
bosonic 50(4) 5 5 states arising from the closed su(2) subsector fall into the —11/3 [280^] level. 
This is the counterpart of the —4 [280#] level in table ^2 To compare with Minahan and 
Zarembo's Bethe ansatz calculation of the corresponding gauge theory operator dimension, 
we must evaluate eqn. (J3.17)) with the appropriate choice of parameters. In particular, 
M n = 2 when two mode indices are allowed to coincide and, comparing with eqn. ([3.29)1 . we 
find perfect agreement with the string theory prediction A = —11/3. States at level L = 4 
in the second multiplet in table El correspond to operators in the sl(2) closed sector of the 
gauge theory and the eigenvalue A = —7/3 [280#] amounts to a prediction for the one-loop 
anomalous dimension of that class of gauge theory operators. As mentioned at the end of the 
previous subsection, we will need to develop a numerical treatment of the sl(2) spin chain 
Hamiltonian in order to assess the agreement between string theory and gauge theory in this 
sector. 

4 Three impurity string spectrum: all orders in A' 

In the previous section, we have studied the eigenvalue spectrum of the string theory per- 
turbation Hamiltonian expanded to leading order in 1/J and to one-loop order in A'. The 
expansion in A' was for convenience only since our expressions for matrix elements are exact 
in this parameter. We should, in principle, be able to obtain results that are exact in A' 
(but still of leading order in 1/J). This is a worthwhile enterprise since recent progress on 
the gauge theory side has made it possible to evaluate selected operator anomalous dimen- 
sions to two- and three-loop order. The simple one-loop calculations of the previous sections 
have given us an overview of how the perturbed string theory eigenvalues are organized into 
superconformal multiplets. This provides a very useful orientation for the more complex 
all-orders calculation, to which we now turn. 

4.1 Inequivalent mode indices: (q ^ r ^ s) 

Our first step is to collect the exact matrix elements of the perturbing Hamiltonian between 
three-impurity states of unequal mode indices. The block structure of the perturbation 
matrix in the spacetime boson sector is given in table H3 and the exact form of the (1,1) 
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block is 



2uJ q LO r U) s 



LUr 



S CD S AF (s 2 + q 2 (l + 2s 2 X')) 

-(q 2 + s 2 )5 cd 5 af - 2qs(5 ad 5 cf - 5 ac 5 df ) + (q 2 - s 2 )5 af 5 c ' d ' - (q 2 - s 2 )5 a ' f '5 cd 
+ + s 2 )5 Cd' 5 a>f> + 2qs (5 a ' d '5 c ' f ' - S a ' c 'S d ' f ')\ + (C^B, D^E, s ^ r) 

+ (A^B, F^E, g^r) j , (4.1) 



where we define u q = ^ q 2 + 1/ A' to simplify this and other similar expressions. 

The off-diagonal H B f contributions to the (1,2), (1,3) and (1,4) blocks are yet more 
complicated. To simplify the expressions, we define 



Fi = yj (u q + q){uj r - r) F 2 = sj (u q - q)(u r + r) 

F 3 = yj (u q - q)(u r - r) F 4 = ^ (u q + q) (u r + r) . 



(4.2) 



Using these functions, the matrix elements in these off-diagonal subsectors are given by: 

5 CD 



(J\a A a B a°(H BF )a^b?W \J) 



32u q u r J I \/~\ 



(F 1 - F 2 )5 AB P 5 - 2(q - r)(F 3 + F 4 )5 AB S^ 



ABryS 



+A(q - r)(F 3 + F^V ~ 2{q + r)(F 3 - F 4 )( 7 ab ')^ 

+ (2qF 3 - 2qF 4 + 2rF 3 - 2rF 4 )( 7 a ' b ) 7<5 - (AqF 3 + AqF A - ArF 3 - ArF 4 )( 1 a ' b ')' 15 

+^=(F 2 - Fi)5 a ' b ' 5 1& + A(q - r)(F 3 + F 4 )5^5 a ' b ' - 2(q - r)(F 3 + F 4 )5^ s (5 ab - 5 a ' b ') 
v A' 

-A\'u q u r (q - r)(F 3 + F 4 )5 AB 5^ + 4v / A 7 (gr - u q u r ) [(F, + F 2 ) {[^'f 5 - ( 7 a ' fe ) 7<5 ) 
-(Fx - F 2 )5^ s (S ab - 5 a ' b ')] + 2(u q + u r )(F 3 + F 4 ) [( 7 a6 ') 75 - W'Y] 



+AV\>(rw q - qtu r )(F 1 + F 2 ) [( 7 afe )^ - ( 7 a ' 6 ')^ 

-AX'(q - r)(F 3 - F 4 ){ru q + quj r )5 AB 5^ - \'5 J 
+2\\u q uj r - qr){q - r){F 3 + F 4 )5 AB 5^ S + AV\'{uj q uj r + gr)(F x - F 2 )6 AB F S 

-2\'(q - r){u q u r + qr)(F 3 + F 4 )5 AB 5A . (4.3) 
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The H FF contribution to the lower-diagonal blocks (2,2), (3,3) and (4,4) is 
(J\b^ r af(R FF )^h^ \J) = 



5 



AB 



48u; r a;o J 



X'< 2rs^l/X' 



((7 ii ) a5 (7 ii ) /37 - (7 i ' j ') a5 (7 i ' j ') /37 j - ^(7 ij ') Q/3 (7^) 7<5 - (ryi'j'yP^'j'y^ 

26 aS 6^(s 2 y/l/X -2rsVXuj r uj s + r 2 (2s 2 VX + V / V^ 7 )) }• (4.4) 



-12 



The bose-fermi Hamiltonian Hbf contributes the following matrix elements to the same 
lower-diagonal blocks: 



(J\b a y r af(H BF )af^ | J) 
+su q {2r 2 \ r X' + 



2u q u r u s 



sVX'5 ab 5 a5 5^ \suj r (2q 2 VX+ ^l/X') 



2u q u r uj s (q + r) v A' 



2uj r q\l + s 2 X') 



+s uj t + 2u! q r (1 + s A') + s uj q — 2s(q + r)X'uj q uj r u s 



+sru q 5 aS 



a'b'\/3-y 



sqtu. 



(4.5) 



To simplify off-diagonal elements in the (2,3), (2,4) and (3,4) blocks, we define 



Gi = a/ (u r + r)(u s - s) 



G 2 = a/ (uj t - r)(uj s + s) 



G 3 = a/ (u r - r)(u s - s) G 4 = A/(cj r + r)(u s + s) . 

The matrix elements in these subsectors are then given by 
(J\b a q b^f(H BF )a^b s ;\J) = 



(4.6) 
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ab ! 



{G l + G 2 ){2-2X'u r u s 



+ (r + s)VX'(g a - A'G 4 (r - u r )(s - u s ) + G 3 (-l + rsX' + ru s X' + u r (s + u 3 )X')j 



+2VA' 



V + s)(G 3 - G 4 ) + VX'(G 1 - G 2 )(rco s - sco r ) ( 7 a T - (l a ' b 'f S 



H-a/A 7 2rsVX'G 1 - 2rsv / A 7 G' 2 + (r - s)(G 3 + G 4 ) + (cj s — u r )(G 3 — G 4 ) 



+2u r uj s VX'(G 2 - GOj [( 7 ab f d - (7 
+ (r + s)v / A 7 ^-G , 4 - A'G 4 (r - u r )(s - u s ) 

+G 3 (1 + rs\' + ru s X' + w r (s + u s )X') 



+ 2^'V 



-2rsv / A 7 (Gi - G 2 



(4.7) 
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The entries in the spacetime fermion block matrix of table |H1 are far too complicated to 
write out explicitly: they are best generated, viewed and manipulated with computer algebra 
techniques. The explicit formulas, along with a collection of the Mathematica programs 
written to generate and work with them, are available on the web |2()| . 

We were not able to symbolically diagonalize the complete perturbation matrix built 
from the exact (in A') matrix elements listed above: with the computing resources available 
to us, the routines for diagonalizing the full 2,048-dimensional matrices would not terminate 
in any reasonable time. As noted in the previous section, however, gauge theory arguments 
suggest that there are three protected SO (4) x 5*0(4) irreps that do not mix with any other 
irreps. It is a straightforward matter to project the perturbation matrix onto these unique 
protected irreps to obtain analytic expressions for the corresponding exact eigenvalues. In 
fact, the superconformal multiplet structure of the three-impurity problem is such that the 
energies/dimensions of all other irreps can be inferred from those of the three protected 
irreps. Hence, this method will give us exact expressions for all the energy levels of the 
three-impurity problem. 

Consider first the sl(2) closed sector. The dual sector is generated on the string theory side 
by bosonic creation operators completely symmetrized (and traceless) on S0(4)Ads vector 
indices. The simplest way to make this projection on eqn. (|4.1|) is to compute diagonal 
elements between the symmetrized states 

a M a? a? \ J) , (4.8) 

with a 7^ b ^ c (and, of course, a,b,c G 1, . . . , 4). The charges of the fermionic oscillators 
under this subgroup are ±1/2, so the three-boson state of this type cannot mix with one 
boson and two fermions (or any other state). Hence, the above projection of eqn. (j4.1|) yields 
the closed sector eigenvalue correction 

5E AdS (q,r,s,J) = - — \qs(l - qsX')u r + qr(l - qr\')u s + rs(l - rsX')uj q 

JLU q Ul r Ul s { 

+ [qr + s(q + r)] X'u q u r LU s ^ 

~ j{-2(g 2 + gr + r 2 )A'-y(A 2 (g + r) 2 )A ,3 + ...| . (4.9) 

To facilitate eventual comparison with gauge theory results, we have performed a small- A' 
expansion in the final line with the substitution s — ► — (q + r) (since the mode indices satisfy 
the constraint s + q + r = 0). The leading correction —2(q 2 + qr + r 2 )X' reproduces the 
one-loop eigenvalue Abb = —2 [280b] located at level L = 4 in the SO(4)Ads multiplet in 
table Ell 

The closed su(2) sector is generated by bosonic creation operators completely sym- 
metrized on traceless SO(4)ss indices. Projection onto this irrep is most simply achieved 
by choosing all mode operators in eqn. ()4.1|) to carry symmetrized, traceless SO (4) 55 labels 
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(they can also be thought of as carrying charge +1 under some SO (2) subgroup of SO (4) 35). 
Direct projection yields the SO(4) S 5 eigenvalue 



5E s5 (q,r,s,J) = -—^ \[qr + r 2 + q 2 (l + r 2 X')] u s + [qs + s 2 + q 2 (l + s 2 \')] 

J UJqUJ r LO s y 

+ [rs + s 2 + r 2 (l + s 2 X')] uj q - [rs + q(r + s)) \'u) q L) r u s j 
« j|-4(g 2 + qr + r 2 )\' + (q 2 + qr + r 2 ) 2 \' 2 

~\ {(f + 3g 5 r + 8gV + llg 3 r 3 + 8q 2 r 4 + 3gr 5 + r 6 ) A' 3 + . . . J . (4.10) 

This is the all-loop formula corresponding to gauge theory operator dimensions in the closed 
su(2) subsector; the leading-order term — 4(q 2 + qr + r 2 ) X' reproduces the one-loop eigenvalue 
Abb = -4 [280b] at level L = 4 in the SO(4) S 5 vector multiplet in table [TT] 

The eigenvalue of the symmetrized pure-fermion irrep can be obtained by evaluating 
the exact matrix element ifpF acting on three symmetrized fermionic creation operators 
with SO (4) x SO (4) indices chosen to lie in the same II projection (with inequivalent mode 
indices). The exact energy shift for this irrep turns out to be 

5E FeTTai (q,r,s,J) = -— — ^ \-4(rs + q(r + s))\'uJ q u r uj s 

4 JuJ q UJ r UJ s { 

(2s 2 + 4r 2 s 2 X' + 2r 2 ) + (s -> r, r -> q, q s) + (q ^ r) 



+ 



(Jjq 



w j|-3(g 2 + qr + r 2 )\' + ^(q 2 + qr + r 2 ) 2 X' 2 

- — (2q 6 + 6g 5 r + 21g 4 r 2 + 32g 3 r 3 + 21g 2 r 4 + 6gr 5 + 2r 6 )A /3 + . . . 1 . (4.11) 
16 J 

The leading-order A' correction — 3(g 2 + gr + r 2 )A' reproduces the Aff = —3 [580^] eigenvalue 
at the L = 4 level in the spinor multiplet in table ^2 This and the higher order terms in 
the eigenvalue will eventually be compared with the dimensions of operators in the closed, 
fermionic su(2|3) sector in the gauge theory. 

The argument we are making relies heavily on the claim that the perturbation matrix is 
block diagonal on the closed subsectors described above: we have evaluated the exact energy 
shift on these subsectors by simply taking the diagonal matrix element of the perturbing 
Hamiltonian in a particular state in each sector. We will now carry out a simple numerical 
test of the claimed block diagonalization of the full perturbing Hamiltonian. The basic idea 
is that, while it is impractical to algebraically diagonalize the full 2, 048 x 2, 048 perturbation 
matrices, it is quite easy to do a numerical diagonalization for a specific choice of A' and 
mode indices q, r, s. One can then check that the numerical eigenvalues match the analytic 
predictions evaluated at the chosen coupling and mode indices. For definiteness, we choose 

q=l r = 2 s = -3 A' = l. (4.12) 
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The predicted eigenvalue shifts of the three protected states, evaluated at the parameter 
choices of (I4.12|) are given in table These values come directly from eqns. (|4. 914. 1014. 11|) 
above (with J set to unity, for convenience). Since we want to compare these energies 



5E : A' = 1 



q = I, r = 2, s = — 3 



8E AdS (l,2,-3,J=l) 
5E s s(l,2,-3, J = 1) 

^Fcrmi(l, 2, —3, J = 1) 



-16.255434067000426 



-20.137332508389193 



18.19638328769481 



Table 15: Exact numerical eigenvalues of three-impurity protected sectors 



to a numerical diagonalization, we must maintain a high level of precision in the numer- 
ical computation. With the parameter choices of ()4.12j) . the numerical diagonalization of 
the full 2, 048 x 2, 048 perturbation matrices on both the spacetime boson (table EJ) and 
spacetime fermion (table EJ) sectors yields the spectrum and multiplicities displayed in ta- 
ble The multiplicities are consistent with the superconformal multiplet structure we 
found in the one- loop analysis (given in table ITT|). The predicted closed sector eigenvalues 
(listed in table IT3j) match, to the precision of the calculation, entries in the list of numer- 
ical eigenvalues. These energies also appear at the expected levels within the multiplets. 
-EUds(l> 2, —3, J) and £"55(1, 2, —3, J) appear in bosonic levels with multiplicity 280#, while 
energy E Fermi (l, 2, — 3, J) appears as a fermionic level with multiplicity 560^; according to 
table E] these are uniquely identified as the central L = 4 levels of their respective multi- 
plets, exactly where the protected energy levels must lie. All of this is clear evidence that the 
'closed sector' states of the string theory do not mix with other states under the perturbing 
Hamiltonian, thus justifying our method of calculating their exact eigenenergies. 

At one loop, we found that the three superconformal multiplets were displaced from each 
other by precisely the internal level spacing. This led to an accidental degeneracy which is 
lifted in the exact dimension formulas we have just derived. To explore this, it is useful to 
have formulas for the eigenvalues of all the levels in each multiplet. ^From the discussion 
in section 3.3, we see that each level in the string energy spectrum can be connected by 
a simple integer shift in the angular momentum J. Since we are working at 0(1/ J) in a 
large- J expansion, all contributions from this shift must come from the BMN limit of the 
theory. In other words, by sending J — > J + 2 — L/2 in the BMN formula for the energy 



we can generate an expansion, to arbitrary order in A', for each level L in the entire super- 
conformal multiplet. 




(4.14) 
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SE(1,2,-3,J= 1) A' = 1 


Mult. 


-30.821354623065 


4b 


-26.9394561816763 


4b 


-26.2093998737015 


64b 


-25.4793435657269 


112b 


-21.5974451243382 


112b 


-20.8673888163637 


448b 


-20.1373325083891 


280b 


-16.2554340670003 


280b 


-15.5253777590258 


448b 


-14.7953214510512 


H2b 


-10.9134230096624 


112b 


-10.1833667016878 


64b 


-9.4533103937133 


4b 


-5.57141195232456 


4b 



6E(1,2, -3, J = 1) A' = 1 


Mult. 


-28.8804054023706 


8 F 


-28.150349094396 


32 f 


-24.2684506530072 


32 f 


-23.5383943450326 


224 f 


-22.808338037058 


224 F 


-18.9264395956693 


224b 


—18.1963832876947 


560b 


-17.4663269797201 


224b 


-13.5844285383314 


224b 


-12.8543722303568 


224b 


-12.1243159223822 


32b 


-8.24241748099347 


32b 


-7.51236117301893 


8b 



(4.13) 



Table 16: All loop numerical spectrum of three-impurity states (q — 1, r = 2, s = —3, A' 
1, J = 1). Left panel: bosons; right panel: fermions 

For the vector SO(A) AdS multiplet, we find 
5E AdS (q,r,J,L) » ±{(L - 6)(g 2 + qr + r 2 )A' - ±(L - 4)(g 2 + gr + r 2 ) 2 A' 2 



3 r 
+ 16 



2(L - 4)g e + 6(L - 4)g 5 r + 5(3L - 14)g 4 r 2 + 20(L - 5)g 3 r 3 

A' 3 



(g 2 + qr + r 2 
16 



+5(3L - 14)g 2 r 4 + 6(L - 4)gr 5 + 2(L - 4)r 6 
5(L - 4)g 6 + 15(L - 4)g 5 r + (50L - 247)gV + (75L - 394)g 3 r 3 
+ (50L - 247)g 2 r 4 + 15(L - 4)gr 5 + 5(L - 4)r 6 ] A' 4 + . . . | (4.15) 



(for convenience in eventual comparison with the gauge theory, the eigenvalues have been 
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expanded to 0(A' 4 )). The corresponding result for the SO(A) s s vector multiplet is 
5E s s(q,r,J,L) » U(L - 8)(q 2 + qr + r 2 )A' - \{L - Q)(q 2 + qr + r 2 ) 2 A' 2 



3 r 
+ 16 



2(L - 6)g 6 + 6(L - 6)g 5 r + (15L - 92)g 4 r 2 + 4(5L - 31)g 3 r 3 



+(15L - 92)<fr 4 + 6(L - 6)gr 5 + 2(L - 6)r fc 



A 



/3 



(q 2 + qr + r 2 ) 



16 



5(L - 6)g e + 15(L - 6)g 5 r + (50L - 309)gV + 3(25L - 156)g 3 r 3 



+ (50L - 309)gV + 15(L - 6)gr 5 + 5(L - 6)r e 



A' 4 + ... 



(4.16) 



Finally, the result for the spinor multiplet is 

5E Fcrrai (q,r,J,L) » - 7)(g 2 + gr + r 2 )A' - \{L - 5)(q 2 + qr + r 2 f\' 2 

2(L - 5)g 6 + 6(L - 5)g 5 r + 3(5L - 27)g 4 r 2 + 4(5L - 28)g 3 r 3 



3 

+ 16 



+3(5L - 27)g 2 r 4 + 6(L - 5)gr 5 + 2(L - 5)r fc 



A' 



/3 



(g 2 + qr + r 2 ) 



16 



5(L - 5)g 6 + 15(L - 5)g 5 r + 2(25L - 139)g 4 r 2 + (75L - 431)gV 



,3„3 



+2(25L - 139)gV + 15(L - 5)gr 5 + 5(L - 5)r e A' 4 + ...}> . (4.17) 



It is important to remember that, to obtain the energies of the states as opposed to the 
energy shifts 5E, we must add the BMN energy of the original degenerate multiplet to the 
above results: 



E 



BMN 



= V 1 + + Vl + A'r 2 + y/1 + X'(q + r) 2 
= 3 + (q 2 + r 2 + qr)X' - \(q 2 + r 2 + gr) 2 A' 2 + 



(4.18) 



We can conclude from the above formulas that all three multiplets have a common internal 
level spacing given by the following function of A' and mode indices: 



A 



,2 



+ -. 



3 r 



16 

(q 2 + qr + r 2 
16 



2q 6 + Qq 5 r + 15g 4 r 2 + 20g 3 r 3 + 15g 2 r 4 + 6gr 5 + 2r 6 



A' £ 



5g e + 15g 5 r + 50gV + 75gV + 50g 2 r 4 + 15gr 5 + 5(r 



A' 4 + ... 



(4.19) 
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We have expanded in powers of A', but an all-orders formula can easily be constructed. The 
multiplets are displaced from one another by shifts that also depend on A' and mode indices. 
We note that the one-loop degeneracy between different multiplets (see table fTTjl is preserved 
to second order in A', but is broken explicitly at three loops. At this order and beyond, each 
multiplet acquires a constant overall (L-independent) shift relative to the other two. 



4.2 Two equal mode indices: (q = r = n, s = — 2n) 

An independent analysis is required when two mode indices are equal (specifically, we choose 
q = r = n, s = — 2n). The all-loop matrix elements are complicated and we will refrain 
from giving explicit expressions for them (though the complete formulas can be found at 
20 ). As in the unequal mode index case, however, exact eigenvalues can easily be extracted 
by projection onto certain protected subsectors. In particular, the energy shift for states 
created by three bosonic mode creation operators with symmetric-traceless SO(4:)Ads vector 
indices (the st(2) sector) turns out to be 

5E AdsM = - J(1 - » a * + 1/A , { fi^l (3 + + M4 + } 



ij-7n 2 A' + n 4 A' 2 -yn 6 A' 3 + ... j 



(4.20) 



The leading order term in the small-A' expansion is the —7/3 [280b] level L = 4 eigenvalue 
in the A 2 multiplet in table EI The energy shift of the SO(4)ss partners of these states 
(belonging to the su(2) closed sector) is 



5E s5 (n,J) = n2X ' -i\ An 2 + (5 + 4n 2 A') + u n (6 + 8n 2 A') \ 

~ l|_ll n 2 A / + g n 4 A /2_ 101^6 y3 + | _ (4 21) 

The one-loop correction corresponds to the —11/3 [280s] level in the Ai submultiplet of 
table EJ As noted above, the protected symmetrized- fermion (su(2|3)) sector does not 
appear when two mode indices are equal. As in the previous section, we can do a numerical 
diagonalization of the full perturbation matrix to verify that the predicted eigenvalues are 
indeed exact and closed, but we will omit the details. 

By invoking the angular momentum shift J — > J + 2 — L/2'm the BMN limit, we can use 
the energy shift of the L = 4 level to recover the exact energy shifts of all other levels in the 
superconformal multiplets of table EI The energy shifts of the vector multiplet containing 
the protected SO(4)Ads bosonic irrep at level L = 4 are given by the expression 

5E AdS (n,J,L) w ^(^(3L-19)n 2 A'-^(9L-38)n 4 A /2 + ^(99L-464)n 6 A /3 
J ^ 2 2 8 

- — (645L - 3160KA' 4 + . . . 1 . (4.22) 
16 
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The shifts of the multiplet containing the protected SO(4) S 5 bosonic irrep are given by 



Once again, we note that in order to get energies, rather than energy shifts, one must append 
the BMN energy of the original degenerate multiplet to these results. Unlike the unequal 
mode index case, there is no accidental degeneracy between superconformal multiplets span- 
ning the three-impurity space, even at one loop in A'. The level spacings within the two 
superconformal multiplets are the same, but the multiplets are offset from each other by an 
L- independent shift (but one that depends on A' and mode indices). 

5 Gauge theory anomalous dimension comparison 

In the previous sections, we have given a complete analysis of the perturbed energy spec- 
trum of three-impurity string states. The 'data' are internally consistent in the sense that 
the perturbed energy levels organize themselves into proper superconformal multiplets of 
the classical nonlinear sigma model governing the string worldsheet dynamics. Since the 
quantization procedure leaves only a subgroup of the full symmetry group as a manifest, lin- 
early realized symmetry, this is by itself a nontrivial check on the consistency of the action 
and quantization procedure. To address the issue of AdS/CFT duality, we must go further 
and compare the string energy spectrum with the anomalous dimensions of gauge theory 
operators dual to the three-impurity string states. 

The task of finding the anomalous dimensions of BMN operators in the limit of large 
i?-charge and dimension D, but finite A = D — R, is greatly simplified by the existence of an 
equivalence between the dilatation operator of M = 4 SYM and the Hamiltonian of a one- 
dimensional spin chain. This correspondence was first proposed by Minahan and Zarembo 
[3], who showed that the anomalous dimension matrix in the 50 (6) subsector of the theory, 
expanded to one loop in the 't Hooft coupling A = g\ M N c , is equivalent to the Hamiltonian 
of an integrable spin chain. Beisert and Staudacher then showed that a more elaborate 
integrable spin chain describes the action of the one-loop anomalous dimension operator 
on the general single-trace operator ^H]- (The complete one-loop dilatation operator was 
derived in |12j.) The key to this development is the fact that the one-loop spin chain 
Hamiltonian has only nearest-neighbor interactions (in the planar large-iV c limit) and is of 
limited complexity. Higher-loop gauge theory physics is encoded in increasingly long-range 
spin chain interactions which generate a rapidly growing number of possible terms in the 
Hamiltonian .5]. Fixing the coefficients of all these terms by comparison with diagrammatic 
computations would be a very impractical approach. Fortunately, Beisert was able to show 
that, at least for BMN operators in the su(2) closed subsector, general requirements (such 
as the existence of a well-defined BMN scaling limit) suffice to fix the form of the spin chain 




J{2 X ' 2 V ' 8 





(4.23) 
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Hamiltonian out to three-loop order [HI • I n this section we will discuss the use of these 
higher-loop spin chains to generate the information we need on the anomalous dimensions 
of three-impurity operators. We will summarize the salient points here, leaving most of the 
details to a separate paper jTTj. 

We have already noted that there are three closed subsectors of BMN operators in which 
impurities taken from a subalgebra of the full superconformal algebra mix only with them- 
selves: we have referred to them as the sl(2), su(2) (both bosonic) and su(2|3) (fermionic) 
sectors. We will focus our attention on these sectors because their spin chain description is 
simple and their anomalous dimensions fix the dimensions of the remaining three-impurity 
operators in the theory. Spin-chain Hamiltonians incorporating higher-loop-order gauge the- 
ory physics have been constructed for the su(2) and su(2|3) sectors but, as far as we know, 
the sl(2) spin chain is known only to one-loop order. 

Although these spin chains are integrable, methods such as the Bethe ansatz technique 
do not immediately yield the desired results for multiple-impurity anomalous dimensions. 
Minahan and Zarembo did use the Bethe ansatz for the one-loop so (6) spin chain (of which 
the exactly closed su(2) system is a subsector) to obtain approximate multi-impurity anoma- 
lous dimensions |4 , but we need results for all sectors and for higher-loop spin chains. As 
mentioned above, the sl(2) spin chain has phenomenological applications and has been ex- 
tensively developed in that context. It is therefore possible that some of the results we need 
can be extracted from the relevant literature. 4 In the end, since we are looking for a unified 
approach that can handle all sectors and any number of loops, we decided that numerical 
methods are, at the moment, the most effective way to extract the information we need 
about gauge theory anomalous dimensions. Since Bethe ansatz equations exist for most of 
the results that are of interest to us, the numerical results obtained here can eventually be 
checked against the Bethe-ansatz methodology: these exercises will be reserved for jllj . 

We begin with a discussion of the bosonic sl(2) sector. For total i?-charge £ (the .R-charge 
is equal to the number of lattice sites £ in this sector), the basis for this system consists of 
single-trace operators of the form 

Tr (V r Z Z*- 1 ) , Tr (V^Z VZ Z^ 2 ) , Tr (V^Z ZVZ Z £ ~ 3 ) , . . . , (5.1) 

where Z is the SO (6) Yang-Mills boson carrying one unit of i?-charge, D is a spacetime 
covariant derivative operator that scales under the chosen sl(2) subgroup of the Lorentz 
group, I is the total impurity number and the full basis contains all possible distributions of 
T> operators among the Z fields. Conservation of various U(l) subgroups of the i?-symmetry 
group ensures that operators of this type mix only among themselves to all orders in the 
gauge theory (as long as we work in the 't Hooft large- N c limit). This gauge theory closed 
subsector corresponds to the symmetric traceless irrep of SO(4)Ads bosons in the string 
theory (states whose energy shifts are given in eqns. (J4.9|) and (J4.2L)jl ). 

The one-loop spin chain Hamiltonian for this sector has been derived by Beisert [T2| in a 
representation where there is a lattice site assigned to each Z field, and each site supports a 

4 We thank A. Belitsky for making us aware of this literature and for helpful discussions on this point. 
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harmonic oscillator whose level of excitation counts the number of D operators acting on that 
Z insertion. The raising operator a\ therefore corresponds to the insertion of a derivative at 
the 2 th lattice site: 



|£> ~ Tr (Z z ) 



[a 



|£) ~ Tr (Zi^TTZZ*-*) 



(5.2) 



The spin chain Hamiltonian is a sum over £ lattice sites of a nearest-neighbor interaction 
which moves excitations between neighbors while keeping the net excitation number (number 
of impurities /) fixed: 



A 



5> 



0[(2) 
k,k+l 



k=l 



)*(4r 



k'=0 



5 k =k' (h(k) + h(n - k)) - 



6 



k^k' 



k-k' 



air 



k' 



(5.3) 

The weighting of different terms by the harmonic numbers h(n) = Y2r=i r_1 * s a general 
feature of spin chains derived from gauge theories and reflects the infrared behavior of gluon 
emission. 5 The impurity number I is the net excitation number of the £ oscillators. The 
interaction propagates impurities around the lattice and assigns special amplitudes to 'col- 
lisions' of multiple excitations on a single site. Our interest is in finding the spectrum of 
energies of three impurities on a large lattice under these dynamics. 

When £ is large and the number of impurities is fixed, it is natural to seek a virial 
expansion for the eigenvalues of H sl ( 2 \ The above Hamiltonian acts on an isolated impurity 
exactly like the usual lattice Laplacian; it can be diagonalized by passing to momentum 
space. The rest of the Hamiltonian amounts to two- and higher-body scattering vertices 
for the single-impurity pseudoparticles, and standard many-body intuition tells us that such 
interaction vertices are suppressed by powers of £ -1 compared to the leading pseudoparticle 
energies. Since we are in fact looking for an expansion of the energy eigenvalues in powers 
of £ _1 , it is useful to rewrite H Bl ^ in terms of momentum space creation and annihilation 
operators a p , a p defined by the usual discrete Fourier transform. Since the definition of H Bl ^ 
is given in terms of its action on states, rather than as an explicit operator, this rewriting of 
the Hamiltonian takes a little work, but the result is simple [22J: 

A 

8T 2 



rxst(2) = ^_ [ 

" 8tt 2 L 



2^4 sin — ala p 
p ~~ 



4E« 



p+q,r+s 



p,q,r,s 



• 2 F . 2 97T . 2 (p + q)n 

sm — sin — + sin 

£ £ £ 



(5.4) 



The pseudoparticle creation and annihilation operators labeled by integer momenta 

p = 0, 1,...,£ — 1 and obey the standard algebra. The ellipses stand for three- and higher- 
body interactions, which are suppressed by even higher powers of £ _1 . For I impurities 



3 We are indebted to A. Belitsky for clarifying discussions on this point. 
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carrying quantized pseudoparticle momenta rij, this leads to an energy formula 



A v ■> . o ^i^T A . , , v 

Sm ~£~ + £3 ^o-body^i, ...,n I ) + ... (5.5) 

To facilitate comparison with string theory, we have reinstated the zeroth-order term for 
the total dimension minus .R-charge of an /-impurity operator. The true eigenvalues differ 
by small corrections from the lattice Laplacian energies of the free pseudoparticles, which 
are labeled by the pseudoparticle momenta n^. These integers are the gauge theory analogs 
of the string mode indices and, to make contact with string theory, we must take the limit 
£ — > oo while keeping the rii fixed. In this limit, the eigenvalues of eqn. (|5.5|) will scale as 

Edim}) = I + -^ (1 ' 2) (R}) + ^£ (1 ' 3) (M) + 0(AiT 4 ) . (5.6) 
This is the leading-order (in A) version of the more general BMN scaling 

E{\J)^E{\/J\l/J)= Eii,2i+j) (4Y J' j , (5-7) 

t=l, j=0 ^ ' 

which, as we have seen in previous sections, naturally characterizes string energy levels. 
In the case at hand, where we are computing energies to 0(A) only, BMN scaling is a 
consequence of the £ -2 scaling of energy eigenvalues that follows automatically from the form 
of the virial Hamiltonian (|5.4jl . Whether the spectra of higher-order spin chain Hamiltonians 
scale with £ in accordance with eqn. ()5.7|) is a very nontrivial question which we will address 
shortly. 

To compare with the corresponding string theory predictions of eqns. ()4.9|) and (|4.20|) . we 
reorganize those results as follows: we reinstate the BMN energy of the degenerate multiplet 
(J4.18|) (expanded to first order in A'); we replace A' with A/ J 2 and replace J by £. This gives 
specific string theory predictions for the large-£ scaling of one-loop anomalous dimensions 
of the AdS, or sl(2), closed sector. As usual, there are two distinct cases: for unequal mode 
indices (q ^ r ^ s = — q — r), we have 

E AdS (q, r, £) = 3 + (£ - 2)(g 2 + r 2 + qr) A + 0(£~ 4 ) , (5.8) 

while for pairwise equal mode indices (n, n, —In) we have 

E AdS (n, £) = 3 + (3£ - 7)n 2 A + (£" 4 ) . (5.9) 

This matches the expected virial scaling of the spin chain eigenvalues displayed in eqn. ()5.6|) . 
with the specific identifications 

= (l 2 + r 2 + Qr) E A ^ = -2(g 2 + r 2 + qr) E^/E^ = -2 (5.10) 
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for q r s = —q — r, or 

ESS = an' B2g = -7»' E'iS/E ( iS = -V3 (5.11) 

for q = r = n and s = —2n. 

To check these predictions, we numerically diagonalize the spin chain for three impurities 
for a sequence of values of £ (up to 30 < £ < 60 in practice). It is a matter of convenience 
whether we construct the £ x £ Hamiltonian matrix using the position space version (|5.3|) or 
the momentum space version ()5.4j) of the Hamiltonian. We then track how the eigenvalues 
evolve as £ varies and fit the data to a general £~ x expansion in order to extract the spin chain 
coefficients E^), E^) (the vanishing of the coefficient E^} i s a check of BMN scaling, 
but is essentially guaranteed here). The results of this exercise, presented in table [T7[ show 
clear agreement between the string and gauge theory predictions. In the fourth column we 
list the ladder of string mode indices that correspond to the succession of eigenstates in the 
gauge theory, according to eqns. (j5.10j) and (|5.11jl . For the low-lying states in the spectrum, 
the numerical agreement is convincing. As one moves up the ladder of energies, higher-order 
l/£ corrections become more important. We would have to generate data up to higher 
lattice sizes and do a more precise scaling fit to improve the agreement at higher levels in 
the spectrum. 

At this point it is appropriate to say a few words about the role of integrability in this 
problem. It was first argued in [23| that the complete GS action of IIB superstring theory 
on AdSz, x S 5 is integrable. Integrability has since taken a central role in studies of the 
AdS/CFT correspondence, as any precise non-perturbative understanding of integrability 
on both sides of the duality would be extremely powerful. Integrability on either side of 
the duality gives rise to an infinite tower of hidden charges that can be loosely classified as 
either local (Abelian) or non-local (non-Abelian) . In the Abelian sector, contact between the 
integrable structures of gauge theory and semiclassical string theory (a subject which was 
first investigated in |2I|) has been made to two loops in A (see, eg., [TUl 123 I2H1 123 )■ (The 
corresponding problem in the non-local sector was addressed to one-loop order in [2*81 I29| .) 
One of the local gauge theory charges, denoted by Q2, can be shown to anticommute in 
the su(2) sector with a parity operator P (to three loops in A), whose action on a single- 
trace state in the gauge theory is to invert the order of all fields within the trace E|. 
Furthermore, Q2 can be shown to connect states of opposite parity. Taken together with 
the conservation of Q2, these facts imply that all eigenstates in the spectrum connected 
by P must be degenerate. These degenerate states are known as parity pairs and their 
existence can be interpreted as a necessary (but not sufficient) condition for integrability. 
The spectrum in table IT7I exhibits such a degeneracy and makes it clear that parity pairs are 
simply distinct states whose lattice momenta (or worldsheet mode indices) are related by an 
overall sign flip. Since the net momentum of allowed states is zero, parity pair states can 
in principle scatter into each other, and their degeneracy is a non-trivial constraint on the 
interactions. As a small caveat, we note that lattice momentum conservation implies that 
mixing of parity-pair states can only occur via connected three-body (or higher) interactions. 
As the virial analysis shows, at the order to which we are working, only two-body interactions 
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E (l,2) 
sl(2) 




7+1,3) 
sl(2) 




,+1,3) / t+1,2) 


String Modes (g, r, s) 


% Error 


1 + 1.2 x 10- 


9 


-2 - 3.1 x 10- 


7 


-2 - 3.1 x 10~ 7 


(1,0, 


-1) 


0.00002% 


3-7.6 x 10- 


9 


-7 + 1.9 x 10- 


6 


-7/3 + 6.3 x 10~ 7 


(1,1, 


-2) 


0.00001% 


3 - 7.6 x 10- 


9 


-7 + 1.9 x 10- 


6 


-7/3 + 6.3 x 10~ 7 


(-!," 


-1,2) 


0.00001% 


4 - 2.8 x 10- 


7 


-8 + 6.9 x 10- 


6 


-2 + 1.7 x 10~ 6 


(2,0, 


-2) 


0.0001% 


7-2.9 x 10- 


7 


-14 + 7.1 x 10' 


-5 


-2 + 1.0 x 10~ 5 


(1,2, 


-3) 


0.0005% 


7-2.9 x 10- 


7 


-14 + 7.1 x 10- 


-5 


-2 + 1.0 x 10~ 5 


(-!,- 


-2,3) 


0.0005% 


9 -4.1 x 10- 


7 


-18 + 1.0 x 10- 


-4 


-2 + 1.0 x 10~ 5 


(3,0, 


-3) 


0.0005% 


12 + 8.4 x 10 


-7 


-28 - 1.5 x 10" 


-4 


-7/3 - 1.2 x 10- 5 


(2,2, 


-4) 


0.0003% 


12 + 8.4 x 10 


-7 


-28 - 1.5 x 10- 


-4 


-7/3 - 1.2 x 10- 5 


(-2,- 


-2,4) 


0.0003% 


13-7.0 x 10 


-6 


-26 + 1.7 x 10- 


-3 


-2 + 1.3 x 10~ 4 


(1,3, 


-4) 


0.01% 


13-7.0 x 10 


-6 


-26 + 1.7 x 10- 


-3 


-2 + 1.3 x 10~ 4 


(-1,- 


-3,4) 


0.01% 


16 - 1.4 x 10 


-6 


-32 + 3.9 x 10- 


-4 


-2 + 2.4 x 10" 5 


(4,0, 


-4) 


0.002% 


19-7.5 x 10 


-6 


-38 + 2.2 x 10- 


-3 


-2 + 1.1 x 10~ 4 


(2,3, 


-5) 


0.01% 


19-7.5 x 10 


-6 


-38 + 2.2 x 10- 


-3 


-2 + 1.1 x 10~ 4 


(-2,- 


-3,5) 


0.01% 


21 - 3.4 x 10 


-6 


-42 + 8.8 x 10- 


-4 


-2 + 4.2 x 10- 5 


(1,4, 


-5) 


0.002% 


21 - 3.4 x 10 


-6 


-42 + 8.8 x 10- 


-4 


-2 + 4.2 x lO" 5 


(-1,- 


-4,5) 


0.002% 



Table 17: Scaling limit of numerical spectrum of three-impurity s[(2) sector at one-loop order 

are present and the parity pair degeneracy is automatic. The same remark applies to the 
string theory analysis to 0(J~ l ) in the curvature expansion. A calculation of the string 
theory spectrum carried out to 0(J~ 2 ) is needed to see whether parity pair degeneracy 
survives string worldsheet interactions; a discussion of this point will be given in |3T)] . 

We now turn to the closed su(2) sector of gauge theory operators, corresponding to the 
symmetric-traceless bosonic 5*0(4)55 sector of the string theory. The operator basis for this 
sector consists of single-trace monomials built out of two complex scalar fields Z and <ft, 
where Z is the complex scalar carrying one unit of charge under the U(l) i?-charge subgroup 
and (f) is one of the two scalars with zero -R-charge, transforming as an 5*0(4) vector in the 
50(6) ~ U{1)r x 50(4) decomposition of the full i?-symmetry group of the gauge theory. 
The collection of operators 

ti^Z*- 1 ), tr(0 / - 1 Z0Z £ - / - 1 ), tr^- 2 ^ 2 ^- 1 ), ... (5.12) 

(and all possible permutations, modulo cyclic equivalence, of the £ factors) forms a basis 
with / impurities and i?-charge equal to £ — I. The anomalous dimension operator simply 
permutes these monomials among themselves in ways that get more elaborate as we go to 
higher loop orders in the gauge theory. An explicit spin-chain Hamiltonian which incor- 
porates gauge theory physics up to three loops has been constructed by enumerating all 
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possible interaction terms and fixing coefficients by demanding BMN scaling behavior of the 
spectrum for large lattice size [HJ E|. We now turn to a numerical analysis of scaling in 
this sector in order to examine the match to string theory predictions at higher loop orders. 

The complete Hamiltonian in this sector will be written as a sum of terms of increasing 
order in the coupling constant A: 

^ su(2) = E(^) n < 2) - ( 5 - 13 ) 

n ^ ' 

The action of the different terms H^ 2 ^ on the operator basis of eqn. ()5.12j) can be built out 
of permutation operators Pij which exchange the fields on the i th and j th sites on a lattice 
of £ sites. Using the compact notation 

L 

{ni, n 2 ,...} = Pk+ni, k+ ni +lPk+n 2 , k+n 2 +l ■ ■ ■ 1 (5-14) 
k=l 

Beisert, Kristjansen and Staudacher jH] find the following explicit forms for the one- and two- 
loop terms in the Hamiltonian (we will return to the question of three- loop terms shortly): 

#2 U(2) = 2({}-{0}) , fff( 2 » = 2(-4{} + 6{0}-({0 ) l} + {l,0})). (5.15) 

Just as in the 5l(2) case, it is clear that for fixed impurity number I and large lattice 
size £, there is a virial expansion of the one-loop energy eigenvalues essentially identical 
to eqn. (|5.5jl . We therefore expect the exact energy eigenvalues to be labeled by integer 
pseudoparticle momenta. 

To obtain the dependence of the spectrum on A, our strategy will be to develop a standard 
Rayleigh-Schrodinger perturbation theory expansion, treating H^ 2 ^ as a zeroth-order Hamil- 
tonian (obtaining its eigenvalues and eigenvectors numerically), then using non-degenerate 
perturbation theory in H s ^ 2 ^ to obtain the next-order corrections (i.e. taking expectation 
values of H^ 2 ^ in the eigenvectors of H% )■ The expansion coefficients of each eigenvalue 
are numbers which depend on the lattice size £ in some non-explicit way: we have to do the 
calculation for many values of £ and perform an extrapolation in £ _1 in order to find the 
information of interest to us but, as we will see, this is not too difficult. The only possible 
obstruction to this program would be a degeneracy in the spectrum of H2 which would 
oblige us to use degenerate perturbation theory. Although the spectrum of is indeed 

degenerate, the higher charge Q2 constrains matrix elements of the perturbing Hamiltonian 
Ht {2) in such a way that a) the first-order perturbation theory calculation can proceed as if 
the spectrum were non-degenerate and b) the degeneracy is maintained to this order. 

These considerations lead us to a numerical scheme involving a series of steps. First, we 
find the eigenvalues and eigenvectors of H^ 2 ^ for the three-impurity Hamiltonian on lattices 
of length up to £ ~ 30 and fit a power series in £ _1 to the eigenvalues. Table IT^l displays 
the values of the coefficients E^?l 1 ^fi ( as defined in eqn. 15. 6|) for the low- lying levels 
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that we infer from this fit. Second, we obtain a series of values for the 0(A 2 ) corrections 
to the eigenvalues by taking the expectation value of the perturbing Hamiltonian H^ 2 ^ 
between the numerical eigenvectors obtained in the previous step. We fit a power series in 
£ _1 to this data to read off the expansion coefficients E^fy for the low-lying levels, with the 
results displayed in table EDI It is important that the 0(A 2 ) data scales as £~ 4 (i.e. that the 
coefficients of lower powers of £~ x vanish to numerical accuracy) as required by BMN scaling 
(see eqn. 1)5.7)1 ) . Since the Hamiltonian was determined in part by requiring this scaling, this 
is perhaps not a surprise: the real test will be the match of the BMN scaling coefficients to 
string theory data. 

To compare with string theory results for the bosonic symmetric-traceless SO (4) gs sector 
eigenvalues, we need to recast eqns. 1)4.10)) and ()4.21)) as expansions in powers of A and £ _1 . 
We denote by E^'" 1 ^ the coefficient of A n £~ m in the large-£ expansion of the string theory 
energies: they can be directly compared with the corresponding quantities extracted from 
the numerical spin chain analysis. The string theory predictions for scaling coefficients, up 
to second order in A, are given in tabled) As usual, the predictions for three-impurity states 



rp(n,m) 

^ S 5 


(q r 7^ s ) (q — r — n) 


£(1,2) 

£(1,3) 

S 5 
P (2,4) 

P (2,5) 


(q 2 + qr + r 2 ) 3n 2 
2(q 2 + qr + r 2 ) 7n 2 
— \{q 2 + qr + r 2 ) 2 — |n 4 
-2(q 2 + qr + r 2 ) 2 -19n 4 



Table 18: String predictions for su(2) scaling coefficients, to two loops 

with unequal mode indices have to be stated separately from those for states with two equal 
mode indices. The level of agreement between string theory predictions and the spin chain 
numerical results is stated in the last column of tables E3 and 1201 The fit is very good for 
low-lying levels and gets worse as we go up in the spectrum. We are confident that, at any 
given level, the fit may be made arbitrarily precise by generating data out to large enough 
lattice size. We take these results as strong evidence that the string theory analysis agrees 
with the gauge theory up to 0(A 2 ) in this sector. 

We now turn to a discussion of gauge theory physics beyond two loops. As it happens, 
the three-loop Hamiltonian can be fixed up to two unknown coefficients («i and a 2 ) by basic 
field theory considerations |2~T) : 

H* o{6) = (60 + 6«i - 56a 2 ) {} + (-104 + 14^ + 96ct 2 ) {0} 

+ (24 + 2«x - 24a 2 ) ({0, 1} + {1, 0}) + (4 + 6«i) {0, 2} 

(-4 + 4a 2 ) ({0, 1, 2} + {2, 1, 0}) - a x ({0, 2, 1} + {1, 0, 2}) . (5.16) 

Originally, these coefficients were determined by demanding proper BMN scaling in the 
theory and that the dynamics be integrable at three loops (by requiring that the Hamiltonian 
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£(1,2) 
su(2) 


E [l,3) 


E (l,3) /E (1,2) 


String Modes (q, r, s) 


% Error 


1 + 1.1 x 10~ 9 


2 - 1.6 x 10~ 7 


2 - 1.7 x 10~ 7 


(1,0, 


-1) 


0.00001% 


3 -9.9 x 1(T 9 


7+ 1.5 x 10- 6 


7/3 + 5.2 x 10- 7 


(1,1, 


-2) 


0.00002% 


3 - 1.2 x 1(T 8 


7+ 1.8 x 10~ 6 


7/3 + 6.2 x 10~ 7 


(-1,- 


-1,2) 


0.00003% 


4-9.7 x 10~ 8 


8 + 1.5 x 1(T 5 


2 + 3.8 x 1(T 6 


(2,0, 


-2) 


0.0002% 


7- 1.6 x 1(T 6 


14 + 2.6 x 10" 4 


2 + 3.8 x lO" 5 


(1,2, 


-3) 


0.002% 


7- 1.7 x 10~ 6 


14 + 2.6 x 1(T 4 


2 + 3.8 x 10~ 5 


(-1,- 


-2,3) 


0.002% 


9 -2.7 x 10~ 6 


18 + 4.9 x 1(T 4 


2 + 5.5 x 1(T 5 


(3,0, 


-3) 


0.003% 


12 + 3.1 x 1(T 3 


27.5 


2.29 


(2,2, 


-4) 


2% 


12 - 1.1 x 10" 2 


29.9 


2.50 


(-2,- 


-2,4) 


7% 



Table 19: Scaling limit of one-loop numerical spectrum of three-impurity su(2) subsector 



7+2,4) 
^u(2) 


7+2,5) 
^u(2) 


p (2,5) / 7+2,4) 


String Modes (q, r, s) 


% Error 


-0.24999997 


-2.000004 


8.00002 


(1,0,-1) 


0.0003% 


-2.24999 


-19.002 


8.445 


(1,1,-2) 


0.01% 


-2.24999 


-19.002 


8.445 


(-1,-1,2) 


0.01% 


-4.00003 


-31.997 


7.9991 


(2,0,-2) 


0.01% 


-12.25 


-97.99 


7.999 


(1,2,-3) 


0.02% 


-12.25 


-97.99 


7.999 


(-1,-2,3) 


0.02% 


-20.25 


-161.8 


7.990 


(3,0,-3) 


0.1% 


-36.10 


-289.6 


8.02 


(2,2,-4) 


5.0% 


-36.18 


-276.4 


7.64 


(-2,-2,4) 


9.5% 



Table 20: Scaling limit of two-loop numerical spectrum of three-impurity su(2) subsector 

H commute with Q 2 up to three loops, and that Q2 anticommute with P); these assumptions 
set CK12 = 0. By studying an su(2|3) spin chain model, Beisert [H] was subsequently able to 
show that independent symmetry arguments, along with BMN scaling, uniquely set a.\ = 
a 2 = (thus proving integrability at three loops). 

The three-loop Hamiltonian Hq 11 ^ can be treated as a second-order correction to H^ 2 \ 
This allows us to numerically evaluate the 0(A 3 ) contribution to the spectrum by using 
second-order Rayleigh-Schrodinger perturbation theory (there is an intermediate state sum 
involved, but since we are doing the calculation numerically, this is not a serious problem). 
There is also the issue of degeneracy but the existence of a higher conserved charge once 
again renders the problem effectively non-degenerate (the details of this argument will be 
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given in a more detailed study of numerical approaches to the spin chain problem 
The resulting three-loop data for large- £ can be fit to a power series in £ -1 to read off the 
expansion coefficients E g ^} 2 y It turns out that, to numerical precision, the coefficients are 
non- vanishing only for n > 5 (as required by BMN scaling). The first two non- vanishing 
coefficients are displayed in table |^ for low-lying levels. 

These results can be compared with string theory predictions derived (in the manner de- 
scribed in previous paragraphs) from eqn. (|4.10|) . and the accuracy of the match is displayed 
in the last column of table EU The important point is that there is substantial disagreement 
with string results at 0(A 3 ) for all energy levels: the low- lying states exhibit a mismatch 
ranging from roughly 18% to 30%, and there is no evidence that this can be repaired by 
taking data on a larger range of lattice sizes. There is apparently a general breakdown of 
the correspondence between string theory and gauge theory anomalous dimensions at three 
loops, despite the precise and impressive agreement at first and second order. This disagree- 
ment was first demonstrated in the two- impurity regime jS], and additional evidence was 
presented more recently in the context of a semiclassical string analysis It is therefore 
perhaps not surprising that the three-loop disagreement is reproduced in the three-impurity 
regime, but it provides us with more information which may help to clarify this puzzling 
phenomenon. 



£(3,7) /£ ,(3,6) 
rj su(2)/ -^311(2) 


String Modes (q, r, s) 


%Error 


16.004 


(1,0,-1) 


33.4% 


14.114 


(1,1,-2) 


18.8% 


14.114 


(-1,-1,2) 


18.8% 


16.037 


(2,0,-2) 


33.6% 


14.272 


(1,2,-3) 


21.7% 


14.272 


(-1,-2,3) 


21.7% 


15.561 


(3,0,-3) 


29.7% 



Table 21: Three-loop numerical spectrum of three-impurity su(2) subsector 

The same exercise can be repeated for the closed su(2|3) fermionic sector, whose string 
theory dual is comprised of pure fermionic states symmetrized in 5*0(4) x 5*0(4) indices 
in either the (1, 2; 1, 2) or (2, 1; 2, 1) irreps (projected onto Il± subspaces). The spin chain 
system is embedded in Beisert's su(2|3) model, where the fermionic sector of the Hamiltonian 
has been recorded up to two-loop order ^j. Since the relevant points of the numerical 
gauge/string comparison have already been made, we will relegate the details of the gauge 
theory side to ^T] and simply state the final results. We also note that Beisert has provided 
us with the fermionic part of the three-loop vertex in the su(2|3) sector. The large- £ spectrum 
of the three-loop contribution will be scrutinized in jTTj, but, based on the existing evidence 
at three-loop order, we do not expect an agreement with string theory. 
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In this sector, the .R-charge and the lattice length are related by J = £ — 1/2. The 
fermionic one- and two-loop string predictions are therefore found from eqn. ()4.11|) to be 



E&l = (q 2 + qr + r 2 ) E^ n = 

E^L = -\{q 2 + qr + r 2 ) 2 fig£i = ~^ + V + ^ • (5-17) 

As noted above, this sector does not admit states with equivalent mode indices. 

The large- £ su(2|3) fermionic spin chain extrapolations are given at one-loop order in 
table The two-loop data are obtained using the same first-order perturbation theory 
treatment described above in the su(2) sector; the results are recorded in table EH1 The 
two-loop spectrum is subject to stronger £ _1 corrections, but the data are still convincing 
and could be improved by running the extrapolation out to larger lattice sizes. The close 
agreement for the low-lying levels corroborates the match between gauge and string theory 
up to two-loop order. 





£(1,2) 
^su(2|3) 




£(1,3) 
sm(2|3) 


£(1,3) /£(1,2) 


String Modes (q, r, s) 


% Error 


H 


- 1.3 x 10 _ 


10 


-1.9 x 10" 8 


-1.9 X 10" 8 


(1,0,-1) 


0.000002% 


4- 


- 1.0 x 10- 


-7 


1.8 x 10~ 5 


4.6 x 10~ 6 


(2,0,-2) 


0.0005% 


7- 


- 2.5 x 10- 


-7 


4.4 x 10~ 5 


6.3 x 10~ 6 


(1,2,-3) 


0.0006% 


7- 


- 2.5 x 10- 


-7 


4.4 x 10~ 5 


6.3 x 10~ 6 


(-1,-2,3) 


0.0006% 


9- 


-3.9 x 10- 


-7 


7.9 x 10~ 5 


8.7 x 10~ 6 


(3,0,-3) 


0.0009% 


13 


-4.0 x 10 


-6 


8.2 x 10~ 4 


6.3 x 10~ 5 


(1,3,-4) 


0.006% 


13 


-4.0 x 10 


-6 


8.2 x 10~ 4 


6.3 x 10~ 5 


(-1,-3,4) 


0.006% 


16 


-2.0 x 10 


-5 


4.1 x 10~ 3 


2.6 x 10~ 4 


(4,0,-4) 


0.003% 


19 


-3.5 x 10 


-5 


7.3 x 10~ 3 


3.8 x 10~ 4 


(2,3,-5) 


0.004% 


19 


-3.5 x 10 


-5 


7.3 x 10-3 


3.8 x 10" 4 


(-2,-3,5) 


0.004% 



Table 22: Scaling limit of one- loop numerical spectrum of three- impurity su(2|3) fermionic 
subsector 
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£(2,4) 


^su(2|3) 


p (2,5) /p(2,4) 
- C/ su(2|3)/- C/ su(2|3) 


String Modes (q, r, s) 


% Error 


-0.25 


-0.99999 


3.99995 


(1,0,-1) 


0.001% 


-4.00006 


-15.990 


3.998 


(2,0,-2) 


0.06% 


-12.251 


-48.899 


3.992 


(1,2,-3) 


0.2% 


-12.251 


-48.899 


3.992 


(-1,-2,3) 


0.2% 


-20.25 


-80.89 


3.995 


(3,0,-3) 


0.1% 


-42.25 


-168.2 


3.98 


(1,3,-4) 


0.5% 


-42.25 


-168.2 


3.98 


(-1,-3,4) 


0.5% 


-64.00 


-254.6 


3.98 


(4,0,-4) 


0.6% 


-90.26 


-359.3 


3.98 


(2,3,-5) 


0.5% 


-90.26 


-359.8 


3.99 


(-2,-3,5) 


0.3% 



Table 23: Scaling limit of two-loop numerical spectrum of three-impurity su(2|3) fermionic 
subsector 
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6 Conclusions 



The BMN/pp-wave mechanism has emerged as a useful proving ground for the postulates 
of the AdS/CFT correspondence. When the full Penrose limit is lifted, a rich landscape 
emerges, even in the two-impurity regime, upon which the string and gauge theory sides of the 
duality have exhibited an intricate and impressive match to two loops in the gauge coupling 
and first nontrivial order in the curvature expansion. While the conditions under which 
agreement is obtained are substantially more demanding in the higher-impurity problem, 
we have shown that this agreement is maintained for three-impurity string states and SYM 
operators. We expect that these conclusions will persist for four or more impurities. Since 
the Bethe ansatz results of Minahan and Zarembo [3] provide an all-impurity prediction in 
the gauge theory, the methods presented here can easily be employed to perform a simple 
check of this statement, for example, in the closed SO (4)^5 sector of the string theory for 
four (or higher) impurity string states. 

Although the two-loop agreement survives at the three- impurity level, we have also con- 
firmed the previously observed mismatch at three loops in the gauge theory coupling. In 
[3*T] it was suggested that this disagreement may be attributed to an order-of-limits problem. 
This notion was made more precise in ^5] , where it was noted that a certain class of long- 
range spin chain interactions may not survive the small-A expansion in the gauge theory. 
A thorough, quantitative understanding of this proposal has yet to be obtained, however. 
In the end, the analyses carried out here will provide an extremely stringent test of any 
proposed solution to this vexing problem. 
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